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ABSTRACT 

A  theory  for  a  method  of  predicting  chorioretinal  burns  from 
nuclear  weapon  explosions  is  presented.  Description  of  a  computer 
program  in  FORTRAN  II  for  use  on  IBM  704,  709,  or  7090  is  given. 

The  calculation  results  in  spatial  and  time  varying  temperature 
profiles  in  the  energy  absorbing  layers  of  the  eye  tissue.  Applying 
tissue  damage  criteria  to  the  results  will  allow  analysis  of  estimates 
of  danger  and  safe  distances  for  viewing  nuclear  explosions. 
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THEORY  ON  RETINAL  BURN  PREDICTIONS 


Making  predictions  on  the  occurrence  of  retinal  burns  in  persons  or 
animals  from  a  nuclear  detonation  essentially  resolves  itself  into  three 
separate,  distinct  problems.  The  first  is  the  production  of  the  thermal 
energy  from  the  fireball  itself j  second  is  transmitting  this  energy 
throu^  the  intervening  medium  of  atmosphere  to  the  eye;  and  third 
is  transferring  this  energy  throu^  the  eye  and  depositing  it  in  the 
layers  of  the  eye  in  -which  heat  and  an  accompanying  tenpcrature  rise  are 
produced  due  to  the  absorption-  of  this  energy.  After  considering  each  of 
these  three  areas  separately,  they  are  put  together  into  one  model  vhich 
integrates  the  particular  -variables  of  each  area.  This  model  attempts  to 
calculate  the  temperatiare  rise  in  the  retinal  layers  for  a  gi-ven  e:q)losion- 
observer  relationship. 

SECTION  I  ■  THEEaiAL  ENERGY  PRODUCTION  ERCM  THE  FIREBALL 

For  he  purposes  of  this  study,  it  is  eussumed  that  the  fireball 
radiates  <s  a  black  body  obeying  Planck's  law  for  the  energy  density  of 
an  isothe  mal  black  body  (Ref.  23): 

f  >  dX  a  0^  .  _ i _ dX  ergs/cm^  (l.l) 

X5  Jch/XkTfi,)_^ 

E  a  l/hcAf  »•  ■  j  Y  ^  ceJ.s/cm^/sec  (1.2) 

where:  k  =  1.380  x  lO"^^  erg/°C  (Boltzmann  Constant) 
h  =  6.625  X  10"^^  erg  sec  (Planck's  Constant) 

X  =  wavelength 
c  a  3  X  10^*^  cm/sec. 
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=  tei^eratm-e  of  radlatlog  tody  (in  this  ease,  the  fireball) 
4*  =  energy  density  (ergs/cm^) 

E  =  total  emissive  power  "between  and  Xg  of  ^  (cals/cm^/sec) 


4.2 


Equation  1.2  states  the  energy  transfer  across  the  surface  maintained 
at  teniperat\are  (T|,]j).  It  is  seen,  therefore,  that  in  order  to  estimate 
the  energy  flux  across  the  fireball,  surface  as  the  explosion  progresses 
in  time,  it  is  necessary  to  estimate  the  temperature  at  i  I’lre- 
ball  is  radiating  from  the  tine  of  explosion  of  the  nacleor  materials  of 
the  bomb  to  the  time  in  which  most  of  the  energy  is  disipated  from  the 
x''ireball.  Jb  should  be  recognized  that  the  phenomena  of  fireball  growth 
depends  upoTi  the  type  of  atmosphere  in  which  it  is  exploded:  if  the  shot 
is  exploded  within  the  atmosphere  the  fireball  behaves  differently  than 
an  explosion  in  outer  space.  Table  1  gives  the  sets  of  values  of 
temperatvires ,  radii,  and  corresponding  times  after  explosion  for  a  typical 
1  KT  explosion  in  sea  level  air  which  have  been  used  in  this  program. 

It  is  assumed  that  the  intensity  of  the  fireball  (temperature)  does  not 
vary  greatly  with  yield,  and  that  the  various  phases  (tine  development) 
of  the  fireball  phenomena  scale  as  the  square  root  of  the  yield.  Also 
given  are  the  corresponding  radii  of  the  fireball  for  the  given  times. 

It  is  assumed  that  these  radii  scale  with  yield  to  the  0.4  power  of  the 
yield. 

In  the  following  paragraphs,  reference  should  be  made  to  Table  2 
which  lists  the  program  of  che  model  in  Fortran  II  language  and  procedure. 
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These  relationships  were  scaled  from  information  in  Chapters  2  and 
7  of  Reference  15 .  The  range  of  time  given  therein  is  not  for  suffl- 
cently  short  or  long  enou^  times;  therefore,  strai^t  line  extrapolation 
on  log-log  plot  has  been  used  to  find  the  hi^er  temperatures  at  shorter 
times,  and  strai^t  line  extrapolation  was  also  used  up  to  two  seconds, 
after  which  a  negligible  temperatxjre  and  a  constant  radius  are  used  to 
fill  out  the  times  (for  computational  purposes).  Because  two  seconds 
is  more  than  60  t^„y,  there  will  be  less  than  lOjt  of  the  thermal 
energy  remaining  in  the  fireball,  and  its  radiation  rate  will  be  very 
slow  compared  to  earlier  times. 

It  is  necessary  to  calculate  the  energy  release  as  a  function  of 

wavelength  as  well  as  time  (throu^  the  temperature  variation)  because 

the  atmosphere  is  wavelength  selective  as  to  transmissivity,  as  is 

also  the  pre-retinal  ocular  media.  In  addition,  recent  research  results 

(Reference  1)  have  shown  the  energy  absorption  in  the  retinal  layers  to 

be  wavelength  dependant,  so  the  energy  must  be  carried  throu^  the 

entire  model  as  a  function  of  wavelength  before  being  absorbed.  As 

a  practical  matter,  the  wavelength  range  is  limited  to  the  range  be- 

o 

tween  3500  and  15,000  Angstroms,  because  of  the  ultra-violet  cut-off 
in  the  atmosphere  and  the  infra-red  limitation  in  the  spectral  absorp¬ 
tion  of  the  pre-retinal  ocular  media.  A  wavelength  band  width  of 
o 

100  A  is  used. 
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There  Is  no  doubt  that  the  interpretation  of  the  fireball  irradiance 
in  this  manner  is  grossly  simplified,  but  because  of  the  level  of  know¬ 
ledge  in  the  other  areas  and  because  of  limitations  of  certain  data,  it 
is  not  deemed  necessary  to  apply  more  sophisticated  techniques  until 
more  complete  information  is  assembled  on  tissue  burn  criteria,  etc. 
However,  modifications  of  the  temperatures  and  radii  used  to  better 
reflect  actual  situations  would  be  well  worthwhile  in  the  prediction 
process  when  such  Information  becomes  available. 

In  the  computation  of  the  energy  emission,  the  prediction  program 

o 

does  it  on  a  spectral  basis  for  each  of  115  wave  bands  of  100  A  width, 
beginning  at  3500 

^  cTx  Ex 

where;  Ix  is  irradiance  on  retina,  uncorrecte'd  for  pre-retinal  ocular 

2 

media  scattering  and  absorption  (cal/cm  /sec/ cm)  (FLL(I,J)  in 
program) . 

c  is  appropriate  units  conversion  factor. 

is  atmospheric  transmission  value. 

2 

Ex  =  fireball  emissive  power  (cal/cm  /sec/ cm) 

Nf  =  "f -number”  of  eye. 

This  information  (the  Ix^^  ELL  (K)  values) together  with  the  fireball 
radius,  is  stored  on  tape,  after  applying  atmospheric  corrections,  for 
the  particular  time  period  under  consideration.  It  should  be  noted  that 
the  program  utilizing  Sense  Switch  #  3,  (Table  2),  will  allow  energies  to 
be  read  directly  into  the  relaxation  program,  so  that  the  method  outlined 
above  for  computing  energy  is  not  a  pre-requisite  to  using  the  program 
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to  compute  the  tissue  temperature  rise.  In  order  to  know  the  energy  per 
unit  area  falling  on  the  lens  of  the  eye  and  within  the  Image  on  the 
retina  (uncorrected  for  ocular  scattering  and  absorption)  the  values  of 
SinO)  and  SUM  are  calculated.  These  values  represent  the  total  energy 
per  unit  area  delivered  to  1)  the  most  center  point  of  the  Image  and  2) 
upon  the  lens  surface,  respectively. 

SECTION  II;  ATMOSPHERIC  TRANSMISSION 

The  amount  of  attenuation  of  radiant  energy  emitted  tr<m  the  fireball 
by  the  atmosphere  depends  upon  the  energy  distribution  (as  a  function  of 
wavelength) ,  the  atmospheric  composition,  and  the  distance  the  observer 
Is  fr(»i  the  source.  For  the  retinal  burn  problem,  any  energy  which  Is 
either  scattered  or  absorbed  does  not  contribute  to  production  of  a 
temperature  rise  In  the  Image  area  or  the  retina,  imless  the  scattered 
radiation  Is  multiply  scattered  In  such  a  manner  as  to  appear  optically 
to  originate  in  the  fireball.  Such  multiple  scattering  Is  neglected  In 
this  model. 

Degredatlon  of  the  energy  Is  mainly  due  to  scattering  of  the  dry 
air  and  scattering  and  absorption  of  the  water  vapor  In  the  air.  Quan~ 
titatlvely,  the  amount  of  scattering  Is  dependant  upon  the  number  of 
scattering  molecules  (air  density)  while  the  amoxmt  of  water  vapor 
determines  the  scattering  and  absorption  (temperatiire  and  relative 
humidity).  Therefore,  an  atmosphere  must  be  chosen  and  the  density, 
temperature,  and  relative  humidity  specified  quantitatively.  The  absorp¬ 
tion  coefficients  of  a  standard  atmosphere  for  each  wavelength  band  must 
be  corrected  for  the  atmosphere  chosen,  then  the  total  transmission 
coefficient  for  each  wave  band  determined  over  the  entire  distance  of  the 
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path  between  fireball  and  the  receiving  eye  (References  17  and  18). 

In  the  part  of  the  prediction  model  dealing  with  this  calculation, 

two  methods  for  calculating  this  transmissivity  are  used,  both  of  which 

are  based  on  a  "flat  earth'*  model.  At  large  distances,  a  "flat  earth" 

model  will  overestimate  the  transmission.  Calculations  of  TRAN  (L)  or 

TRAN  (M)  (Table  3)  depends  upon  whether  there  is  an  altitude  difference 

between  the  source  and  receiver.  If  an  altitude  difference  does  exist, 

TRAN  (L)  is  used:  if  there  is  not,  TRAN  (M)  is  calculated.  The  basic 

difference  is  that  TRAN  (L)  allows  for  a  variation  in  atmospheric 

density  in  the  choice  of  constants  and  in  the  method  of  calculation 

(based  on  concept  of  "reduced  height"  atmosphere) .  These  transmission 

coefficients  are  calculated  for  each  of  the  115  wavelength  bands  of 
o 

100  A  width  and  stored  in  memory  to  be  used  to  calculate  the  energy 
emission  spectrum  (FliL(K)),  described  in  Section  I. 

In  dealing  with  the  calculation  of  the  atmospheric  transmission, 
the  situation  that  develops  may  be  described  as  being  "extremely  precise 
and  highly  Inaccurate".  An  intricate  calculation  method  i  been 
developed  and  is  presented  for  use.  However,  because  of  the  limited 
accuracy  of  the  input  data  and  because  of  the  fallacy  of  trying  to 
estimate  temperatiire ,  relative  humidity,  and  atmospheric  composition 
under  an  average  situation,  the  power  of  the  method  is  limited  to  the 
accuracy  of  the  known  and  specified  conditions.  Indeed,  here  is  one 
area  where  a  general  method  cannot  overcome  lack  of  specific  data  and 
facts. 

The  following  paragraphs  will  outline  the  method  developed  and  data 
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will  be  given,  such  as  It  is,  for  a  typical  Pacific  Ocean  atmosphere  on 


a  typical  winter  day  at  a  typical  time.  It  should  be  stressed  that  lack 
of  experience  with  the  program  precludes  a  knowledgeable  evaluation  of 
the  effect  of  various  changes  in  the  atmospheric  input,  but  it  is  felt 
that  unless  appreciable  range  of  values  for  a  specific  location  or 
situation  exist,  the  atmospheric  transmission  may  be  regarded  as  a 
second  order  effect  in  view  of  the  highly  Influenclal  uncertainties  in 
the  other  areas,  such  as  weapon  radiation  output  and  tissue  temperature 
damage  criteria. 

Method 

The  method  used  to  estimate  the  atmospheric  transmission  is  that 
outlined  on  page  432  of  Reference  IS.  The  method  is  based  on  observations 
of  the  transmission  of  solar  radiation  through  a  cloudless  atmosphere 
(assumed  dust  free) : 

^  “  ^aX 

where:  a.y  is  the  spectral  transmission  coefficient  for  a  vertical  path 
(optical  air  mass  imity,  when  P=Pq)« 

a^^j^  is  the  vertical  transmission  considering  only  the  effect  of 
scattering  by  pure  dry  air. 

a^^  is  the  transmission  considering  only  the  effect  of  scattering 
by  1  cm  of  precipitable  water  vapor. 

w  is  the  amount  of  precipitable  water  in  the  path,  measvired  in  cm. 

For  this  study  the  relationship  was  modified  to  accotint  for  any 
straight  path: 

”  ("ax)"  (■Vx)'" 
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where:  is  the  toteO.  transmission  factor  for  any  chosen  c^tical  path, 

m  is  the  number  of  air  masses  in  the  optical  path. 
w*m  is  the  total  precipitable  water  vapor  in  the  chosen  optical 
path. 

Table  3  gives  values  of  a^^j^  and  a^^^  (called  PTRAN  and  QTRAN  in  the 
program)  based  on  Reference  l8,  pages  431>  f.f. 

The  Smithsonian  tables  also  recognize  absorption  by  the  water  vapor 
in  the  air  and  present  corrections  for  this  based  on  a  spectrally  inte¬ 
grated  transmission  coefficient.  This  correction  is  not  applied  to  this 
problem  for  the  following  reasons: 

1.  It  is  not  given  spectrally. 

2.  The  largest  effect  of  the  water  vapor  absorption  would  be 
in  the  lower  and  upper  wavelength  regions  of  the  spectrum  considered  in 
this  study.  Because  these  regions  are  highly  affected  by  the  atmos¬ 
pheric  scattering  and  the  pre-retlnal  ocular  absorption  and  scattering, 
its  total  effect  on  the  total  energy  transmitted  would  be  almost 
negligible,  when  these  other  attenuations  are  considered  and  accounted. 


the  distance  through  an  amoxmt  of  sea  level  atmosphere  equivalent  to  a 

verticsQ.  path  (optical  air  mass  unity,  when  p  =  p^). 

This  is  written  in  the  program  as: 

r  -]s  r  nioeji  5 

TRAN(M)  «  PTRAIJ(M)  QTRAN(M) 
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The  above  equation  has  been  -written  in  the  program  as ; 

TRAK(L)  a  (1.0  -  CON  l^(l.0  -  PTRAN(L)))^  (QTRAN(L))^°®^  ^ 

■where:  A  =  1  + 

D  is  the  horizontal  distance. 

h  is  the  vertical  distance  be-tween  explosion  and  observer, 

represents  the  average  precipitable  water  -vapor  o-ver 
the  slant  distance  per  horizontal  mile  of  optical  path  (measxuced 
in  centimeters  of  water). 

Wjjp  is  the  water  vapor  multiplication  factor  -which  accounts  for 
the  assvmptlon  that  there  is  negligible  water  vapor  above  a  level 
of  5  miles  altitude  above  sea  level  and  that  the  lower  of  the 
observer  or  explosion  altitude  may  or  may  not  be  at  sea  level. 

For  calculation  of  this  factor,  see  Appendix  A. 

C  is  a  constant  which  represents  the  centimeters  of  precipitable 
water  vapor  per  horizonteJ.  mile  of  air  at  sea  level.  For  Pacific 
atmosphere  at  25°  C  and  lOOjt  relati-ve  htmddity,  an  approximation 
for  C  at  sea  le-vel  is  3.W  cm.  (See  Table  in  Appendix  A  for 
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otfrer  values  of  C,  as  a  function  of  altitude.  If  the  type  of 
atmosphere  is  changed,  other  values  of  C  for  the  new  location 
will  he  needed. ) 

is  the  relative  humidity  averaged  with  water  content  over  the 
cqptical  path. 

D  is  the  horizontal  distance  between  observer  and  explosion. 
SECTION  III:  ENERGY  TRANSMISSION  AND  ABSORPTION  IN  THE  EXE,  AND  METHOD 
OP  CALCULATING  THE  THiPERATURE  PRODUCED 
(The  calculation  of  a  tenperat\nre  rise  in  the  retinal  layers) 
Research  results  (Reference  l)  have  yielded  values  which  may  be  used 
as  energy  transmission  values  (as  a  function  of  wavelength)  throu^  the 
pre-retlnal  ocular  media.  These  are  presented  as  a  function  of  the  11^ 
wavelength  bands  of  100  ^  each  in  Table  3>  Table  ^  presents  the  total 
absorption  coefficients  (in  decimal  percent)  used  for  the  layers  of 

plgnented  eplthellxim  and  the  choroid,  as  a  function  of  wave  length  for 

o  0 

100  A  bandwldths,  starting  at  33^00  A  (Reference  l).  !lhese  three  factors 


are  appropriately  applied  to  determine  the  absorption  in  each  of  the  two 
layers.  (Refer  to  Figure  1:  1  denotes  region  of  P.E.,  2  denotes  region 


where;  RET  is  thickness  of  P.E.  («  lOp). 
A  X  =  100  A. 
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Absorption  In  Choroid  (Region  l)  —  cals/cau  /sec 

15,000  1 

Q  (l-ABSMO(X))  (I-ABSI(X))  (ABSII(X))  A  X 

(3.2) 


X  «  3500  A 


where  CHOR  is  the  thiclcness  of  choroid  (for  htanans  a60fx  to  lOOp). 


A  correction  factor  of  _ _  — . —  is  applied  to  AI  and  All  if  the 

(3.5  X  10-^) 

radius  of  the  image  (R  in  centimeters)  is  less  that  3.5h,  which  is  considered 
to  be  an  average  image  limitation  value,  based  on  chromatic  abberation,  dif¬ 
fraction,  and  corrected  for  a  — i —  energy  distribution  (Reference  19) . 

X 

The  resulting  values  of  Al  and  All  are  then  used  as  heat  sources  In  a 
heat  conduction  problem  which  attempts  to  calculate  the  spatial  and  time 
varying  temperature  profiles  in  the  volume,  Including,  and  surrounding  the 
volume  defined  by  the  fireball  optical  image  on  the  retina  and  the  layer 
thickness  of  the  absorbing  tissues. 

Certain  basis  assumptions  are  made  about  the  "heat  generation"  model.  It 
is  assumed  that  all  tissue  thermal  properties  are  nearly  the  same  as  water  and 
that  absorption  and  conduction  occur  imiformly.  No  corrections  are  made  for 
changes  in  tissue  temperature.  All  heat  that  is  not  absorbed  in  the  pigmented 
eplthelliun  or  choroid  is  discarded  as  not  causing  a  sigpiif leant  contribution 
to  the  temperature  rise  in  those  two  regions.  There  is  assumed  conduction  in 
the  tissue  surrounding  the  absorbing  voltune,  but  no  heat  generation.  The 
absorbing  volumes  "grow"  in  radius  as  the  image  radius  expands  with  time  as 
a  result  of  fireball  growth,  but  the  thicknesses  of  absorbing  properties  do 
not  vary  with  time  or  change  with  temperature.  Be-radiation  of  energy  due  to 
the  local  temperature  rise  in  the  area  of  the  absorbing  volume  is  neglected. 
There  is  some  indication  (Reference  20)  that  the  absorbing  of  the  radiant 
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energy  In  the  pigmented  epithelium  does  not  take  place  uniformly,  but  Is 
absorbed  in  pigmented  cells  which  are  interspaced  between  semi-transparent 
tissue.  This  local  absorption  effect  is  neglected  in  the  model,  as  quanti¬ 
tative  information  on  this  method  of  absorption  is  not  available.  However,, 
the  .surface  to  voliune  ratio  of  these  cells  is  considered  to  be  so  large  that 
there  should  not  be  appreciable  build-up  in  temperature  within  them  before 
conduction  will  dissipate  the  heat  to  the  surrounding  tissues.  Should  the  heat 
generation  perform  more  as  surface  heating  than  volume  heating  in  these  cells, 
an  entirely  different  method  of  treating  the  heat  source  would  be  required. 
However,  this  remains  to  be  proven  experimentally  or  otherwise  demonstrated 
formally. 

The  fundamental  heat  conduction  equation  is  (Reference  21): 

^  ^  "'ft')  (3.3) 

k  k  St 


vheres  k  =  conductivity  of  material. 

p  s  density  of  material. 

Cp  =  specific  heat  of  material. 

2 

V  =  Lapl6u:lan  operator, 
v^^x  =  temperature  as  a  ftmctlon  of  time, 
t  =  time. 

=  volume  absorption  rate  in  region  (i). 


Equation  3.3  is  put  into  difference  equation  form  according  to  the 


following  diagram: 


t 


J 


I 
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Expressing  the  Laplacian  operator  in  terms  of  a  cylindrical  differ¬ 
ence  equation  for  the  tenperature  of  point  1  in  terms  of  points  adjacent 
to  points  3>  5,  6,  and  9  yields: 

2  2  2 

V  'Vss  o  V  +  1  Sv  +  V  for  S'v  -  n 


^  t  V,  1.  Y3  -T;  4-  Vit  -  2v  j.  Yg 


k 


(  for  region  1) , 


(3.4) 

(3.5) 


**  _ (si  “  (for  region  2), 

k  k 

~  — ^  (for  region  3) 


(3.5) 

(3.5) 

(3.6) 


where:  subscript  n-1  refers  to  the  temperature  at  point  1  at  the 
preceeding  relaxation  time. 
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Equation  3 A,  3.5>  and  3*6  then  are  -used  to  approximate  Equation  3.3  in 


finite  difference  form: 


v^-  2v^  + 


V  - 

3 _ 


V.  -  2v  +  V 


ts?- 


2  +  PI 


V  -  V 

_ j-(° 

*  -  ’'(n  -  l) 


from  which 


V  +  -  V 

_5 _ 3_  +  _A 


Af.^ 


2rZr 


5  ’'2 

+  - 


+  pi  +  £2e  1  -  1) 


t  -  t(n  -  1) 


(3.7) 


/  _ ,2-^2 _ \ 

k(ti  -  t^jj  _  2^))  ^2  ^2  j 


InitieJ.  and  bo\jndary  conditions 

At  t  =  0,  all  grid  points  v(l,j)  «  0, 
where;  1  <  J  <  LAT,  1  <  I  <  LONG,  I  and  J  integers,  LAT  and  LONG 
defining  the  nvmiber  of  grid  points  in  the  radial  and  axial 
directions . 

Equation  3.7;  therefore,  calculates  the  teaa>erature  rise  above  the 
ambient  body  teiig>erature ,  which  may  be  chosen  to  reflect  the  body  teiiq>- 
erature  of  subject  being  considered.  (However,  v(l,J)  ■  T  where  T  is 
any  value,  by  changing  statement  T  ■  0.0  following  statement  11  of  the 
program;  calculated  temperatures  wovild  then  be  actuEJ.  temperatures  if 
T  were  body  temperature . ) 

Absorption  function  PI  is  defined  by  eqtiatlons  (3.5). 

At  the  bo\indaries  between  regions  1  and  2,  1  and  3;  and  2  and  3; 


'2^  ''l 


Ik 


Also ; 


Sv  av^ 


av3 

Sr 

ar' 

az 

az 

^^3 

a^ 

^3 

Sr 

ar  ^ 

az 

bz 

at  the  appropriate  houndries. 


In  region  3,  the  outer  cylindrical  boundary  is  assumed  to  be  non- 
conductive  so  that: 


at  outer  grid  points. 

To  effect  this  condition, 

v(l,J)  -  v(2,J)  ~1 

(t  fr  ok^  L  (Numbers  24,  25,  48,  49  would  depend  upon 

v(I,25;  ■  v(I,24;  ^  __  la^^  Lqjjq 

V(lt9,j)  .  v(48,J)J  ■ 

is  dictated  by  the  program.  The  physical  effect  of  this  is  that  no 
heat  is  lost  from  region  3.  As  long  as  region  3  is  kept  large  coo5)ared 
to  regions  1  and  2,  it  is  felt  that  no  loss  of  heat  from  region  3  coun¬ 
ter-balances  neglecting  absorption  in  region  3*  Heat  generation  in 
the  three  regions  is  governed  by  equation  3.5  stated  previously. 

From  the  above,  it  is  seen  that  the  temperature  of  any  i>oint  in  a 
grid  can  be  found  if  the  heat  source,  the  previous  temperature  of  that 
point,  and  the  surrounding  foiar  temperatures  are  knovm.  A  grid  is  set 
up  in  the  retinal  layers  to  represent  a  radial  plane  of  a  ri^t  circular 
cylinder. 
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II 

CHOROID 


80  -  lOOii 


Figure  1.  Schematic  Diagram  of  Absorption  and  Conduction  Volume 
It  is  assumed  that  the  energy  absorption  (and  hence  the  temperatures) 
does  not  vary  with  angle  ♦ . 

The  resulting  grid  is  as  follows  (numbers  refer  to  sequence  of 
applying  equation  (3.7))» 


9  10  '  11  Is 

§  5  ‘V  5  ih 


1  2  15 


I 

Figure  2.  Relaxation  Grid  Numbered  to  Show  Progression  of  Applying 
Equation  (3.T)« 


Equation  3.7  is  applied  first  at  point  1,  then  at  point  2,  etc., 
in  a  hemispherical  motion  about  a  center  of  point  1.  By  progressing 
in  this  maimer  throu^  the  grid  several  times,  the  temperature  of  each 
grid  point  "relaxes"  to  values  which  approximate  the  temperatures  for 
the  time  specified.  A  test  is  built  into  the  program  so  that  when  the 
largest  change  in  temperature  of  any  point  in  the  grid  is  less  than  a 
pre-assigned  percentage  of  the  largest  temperature  in  the  grid,  the 
relaxation  procedure  is  halted.  The  smaller  this  TEST  percentage  is 
made,  the  greater  the  accuracy  of  the  resxilts.  Experimentation  with 
this  value  is  needed  in  most  situations.  The  grid  in  Figure  2  is 
superimposed  on  Figure  1  so  that  the  starting  i>oint  (point  l)  is  the 
center  of  the  pigmented  epithelium  layer,  on  the  z  axis.  The  size  of 
the  grid  is  chosen  to  be  about  3  times  the  largest  image  radius  and 
about  three  or  four  times  as  long  as  the  thickness  of  the  choroid. 

After  a  relaxation  is  completed  for  one  time  period,  time  is 


incremented  and  new  values  are  calculated  and  the  relaxation 
process  commenced  again.  i5ie  program  halts  or  goes  to  the  next 
program  when  zeroes  appear  in  FTIMS(l)  and  DELTl(l).  (Therefore, 
FTIMS(9)  and  DELTl(9)  must  always  be  zero.) 

Typical  variable  data  input  is  listed  in  Table  6.  Results  of 


the  calculations  for  the  tenperature  profiles  at  various  stages  of  the 
explosion  development  are  listed  in  Table  6. 


Grid  size 

The  grid  tised  presently  for  this  program  is  points  in  the  z 
(axial)  direction  and  25  points  in  the  r  (radial)  direction  (a  point  is 
denoted  by  (l,J)  variables  sis  to  point  location,  i.e.,  (25, l)  denotes 
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center  of  pigmented  epithelium  on  z  axis,  (4,13)  denotes  the  4th  radial 
point  from  the  center  line  vhich  is  in  z  coliamn  number  13. 


The  physical  side  of  the  grid  depends  upon  the  values  for  and  These 
values  are  usueG-ly  selected  by  estimating  the  largest  image  radius,  multi¬ 
plying  this  value  by  3  to  allow  for  conduction,  and  then  dividing  it  by  the 
number  of  radial  space  Increments. 

3^i  (largest)  (largest) 

^  =  8 

It  is  generally  convienient  to  let  ^  but  the  product  of  £sz 

and  the  number  of  axial  grid  points  must  generally  be  at  least  300  to  400 
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microns  in  order  to  eillaw  proper  conduction  axially. 

Provisions  for  other  attenuators 

A  Constant  Attenuation  Factor  (abbreviated  CAF  in  the  program)  is 
provided  so  that  attenuators  such  as  plexiglswss  or  fixed,  constant 
density  (with  wavelength)  filters  may  be  accommodated.  The  proper 
transmission  value  of  the  attenuator  (in  decimal  percent)  is  then  stated 
in  the  input  variable  for  CAF.  If  no  such  attenuator  is  present  in  the 
physical  problem,  a  veilue  of  CAF  »  1.0  is  Inserted  as  the  value  for  CAP. 
(Actually,  the  words  "constant  attenuation  factor"  are  a  misnomer  —  it 
should  be  "constant  transmission  factor"  because  of  its  method  of  use.) 

Calculation  of  Henrlques*  Damage  Integral 

There  is  a  provision  in  the  program  (c.f.  statements  150  to  15^) 
for  calculating  the  integrated  cellular  damage  based  on  the  Henriques* 
Damage  Integral  eniperical  relationship.  This  calculation  can  be  utilized 
by  placing  Serjse  Switch  4  in  the  down  position.  The  calculation  for  H 
is  made  for  variovis  grid  points  which  are  determined  by  reading,  as  input, 
various  values  of  IiEX(l);  a  pair  of  LEX(l)  values  corresponding  to  a 
grid  point  in  terms  of  (l,J).  By  changing  the  statement  following  state¬ 
ment  400  (Signia(L)  »  etc.)  any  relationship  for  Henriques'  function 
may  be  utilized.  It  is  anticipated  that  this  portion  of  the  program 
may  become  increasingly  valuable  in  applying  damage  criteria  as  the 
program  is  utilized  in  analyzing  laboratory  produced  permanent  damage 
where  the  damage  producing  temperature  may  be  estimated  by  staining  tech¬ 
nique.  Ultimately,  however,  uhe  dependence  will  be  directly  related  to  the 
confidence  in  the  Henriques'  relationship  ^^sed. 
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Damage  criteria 

In  or^er  for  this  model  to  he  effective,  some  criteria  must  he 
established  to  determine  when  irreparable  damage  occurs.  Previous 
attempts  (References  22  and  l6)  have  been  to  specify  the  damage  thres¬ 
holds  in  terms  of  irradiance  (cals/cm^/see)  radiant  exposure  (cals/cm^), 
image  size,  and  time  of  exposure.  It  is  now  thou^t  that  the  temperature 
of  the  absorbing  and  conducting  tissue  of  the  retina  is  more  fundamental 
to  the  damage  mechanism  that  the  total  energy  or  energy  rate  of  deposition. 
If  this  be  tme,  the  damage  criteria  should  be  related  to  the  temperature¬ 
time  history  of  the  irradiated  material.  Research  efforts  (at  the  Medical 
College  of  Virginia)  are  now  in  progress  to  attempt  to  demonstrate  this 
point.  If  this  damage  relationship  can  be  fouiid  experimentally,  it  would 
provide  damage  criteria  for  the  predictions  produced  by  this  program. 

Until  the  laboratory  results  are  available,  an  interim  method  to 
estimate  the  bum-damage  criteria  has  been  used  to  evaluate  the  results 
of  calculations  based  on  this  model.  It  is  reasoned  that  irreparable 
damage  occurs  where  the  tissue  is  heated  above  the  tenperatxire  at  which 
boiling  of  the  fluids  will  occur  (  ~  1CX)°  C  or  a  tejcoperature  rise  of 
~  65°  C).  It  is  also  thought  that  irreparable  damage  can  occur  when  there 
is  no  evidence  of  boiling  or  vapor  eruption  in  the  tissue,  which  means 
a  macroscopic  maximum  temperature  somewhat  below  100°  C  was  obtained. 

On  the  lower  end  of  the  scale,  it  is  thought  that  irreparable  damage  is 
not  likely  if  a  temperature  rise  of  "less  than  20°  C  exists  for  a  short 
time  (less  than  a  second). 

o 

Thus,  it  is  reeisoned  that  a  temperature  rise  of  greater  than  20  C, 

.  o 

but  less  than  65  C,  for  a  time  period  of  the  order  of  the  blink  reflex 
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or  less,  will  be  the  tenIperat^lre  range  in  which  irreparable  damage  is 
produced.  A  more  specific  criteria  awaits  further  research  results. 

Results 

Input  values  and  results  and  print-out  for  a  typical  computer  run 
are  included  in  Table  6. 

Operation  Notes 

This  program  has  been  operated  on  IBM  704,  709,  and  7090  computers 
with  32,000  word  memory  units.  Running  time  for  an  average  problem  is 
approximately  20  -  30  minutes  on  the  704  and  5-10  minutes  on  the  7090. 
A  calculation  using  one  yield  and  five  distances  to  determine  a  safe 
distance  requires  about  35  minutes.  Of  course  these  times  are  highly 
dependant  on  the  grid  size,  time  increments,  and  TEST  values  used.  Use 
of  Sense  Switch  #  2  allows  consecutive  runs  to  be  made  without  stopping 
the  computer;  however,  the  same  input  data  must  be  used,  and  only  the 
variable  cards  can  be  changed  from  run  to  run  in  this  procedure. 
Computer  luilts  with  8k  memories  may  be  used  by  utilizing  a  smaller 
number  of  grid  points. 
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APPENDIX  A 


Calculation  of  W  .  (WMFl) 
mr 


Ceilculatlon  of  the  precipitahle  vater  vapor  multiplication  factor 
may  he  demonstrated  from  the  following; 


TABLE  A-1 


Values  of  Precipitahle  Water  Vapor  for  a  Typical  Pacific 


Atmosphere  as  Might  he  Found  Near  Hawaii  is. 

the 

Altitude 

Tenip 

(°c5 

Water 

Relative 

Cm.  of 

Cm  of 

(mi) 

Vapor 

Density 

Water 

Humidity 

Precipitahle 
Water  Vapor 
per 

Precipitahle 
Water  Vapor 
per 

Altitude 

Temp 

(^C) 

Vapor 

Relative 

Horizontal 

Horizontal 

(mi) 

Density 

Humidity 

Mile 

(100^  P.el. 
Humidity) 
Cbnstant  C 

Mile 

at  Indicated 
Relative 
Humidity 

0.0 

Ji5 

0.02178 

0.80 

3.48 

2.78 

17 

0.01447 

U.85 

2.32 

1.98 

1 

Xe 

13 

0.011.35 

0.80 

1.82 

1.46 

10 

0.00941 

0.55 

1.51 

0.83 

2.0 

7.5 

O.OISD*! 

0.42 

1.28 

0.54 

2.5 

4 

0.00636 

0.42 

1.00 

0.42 

3.0 

o.c 

0.00485 

0.42 

0.78 

0.33 

3.5 

-  6 

O.CO35O 

0.42 

0.56 

0.24 

4.0 

-13 

0.00250 

0.42 

o.4o 

0.17 

4.5 

-23 

0,00170 

0.42 

0.27 

0.11 

5.0 

-34 

0.00120 

0.42 

0.19 

00  0 

•  « 

Average  from  0  to  5  miles  ss 


=  0.886  cm  water/mile. 


(A.l) 


0.886 

0*% 


(A.2) 


for  values  of  C  =  3.W  and  H^.  =  0.5, 

=  _ 0-^ _ =  0.32  (A-3) 

3.W  X  0.8 

This  Is  the  value  for  for  a  slant  path  between  altitudes  of 
0  and  5  miles,  and  an  average  relative  humidity  of  805^  in  the  first  several 
thousand  feet  of  atmosphere. 

Should  the  higher  altitude  be  above  5  miles,  a  correction  factor 
must  be  applied  to  to  account  for  the  fact  that  the  entire  horizontal 

distance  (D)  does  not  apply: 


(A.l^) 


If  the  lower  altitude  is  above  sea  level,  the  proportion  of  the 

last  column  in  Table  A>1  should  be  omitted  in  computing  the  average 

precipitable  water  vapor.  For  example,  should  the  lower  altitude  be 

3000'  and  the  hi^er  altitude  5  miles,  the  avereige  precipitable  water 

vapor  would  be  8.86  -  2.78  “  0.67  cm  water/mile.  Kierefore,  equation  A.l 

9 

should  now  be: 

C.Wjnf. %  »  0.67  (A.5) 
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TABLE  1 


.  8 

Normalized  time  vs.  fireball  temperature  and  radius.  (XNORT  x  10  seconds  — 
TEMPI  X  10”  degrees  C  —  RADII  x  10  feet) . 


00000010 

00000020 

00000030 

00000040 

00000050 

00000060 

00000070 

20000 

10000 

9000 

8200 

7600 

7000 

6600 

12 

15 

16 

19 

22 

24 

26 

00000080 

00000090 

00000100 

00000200 

00000300 

00000400 

00000500 

6400 

5800 

5600 

3700 

2800 

2000 

1700 

27 

29 

30 

42 

49 

59 

65 

00000600 

00000700 

00000800 

00000900 

00001000 

00002000 

00003000 

1300 

1100 

920 

800 

740 

340 

210 

69 

75 

80 

87 

89 

120 

140 

00004000 

00005000 

00006000 

00007000 

00008000 

00009000 

00010000 

150 

120 

92 

80 

62 

60 

54 

170 

190 

200 

220 

240 

250 

260 

00020000 

00030000 

00040000 

00050000 

00060000 

00070000 

00080000 

23 

15 

11 

9.0 

7.4 

6.4 

5.5 

300 

430 

480 

530 

560 

580 

620 

00090000 

00100000 

00200000 

00300000 

00400000 

00500000 

00600000 

4.8 

4.4 

2.2 

1.8 

1.9 

2.2 

2.5 

660 

680 

880 

990 

1080 

1140 

1210 

00700000 

00800000 

00900000 

01000000 

02000000 

03000000 

04000000 

2.7 

3.2 

3.5 

3.8 

6.6 

8.0 

8.5 

1260 

1320 

1350 

1360 

1680 

1860 

1960 

05000000 

06000000 

07000000 

08000000 

09000000 

10000000 

20000000 

8.5 

8.2 

7.8 

7.0 

6.4 

5..9 

3.4 

1990 

2040 

2080 

2100 

2160 

2190 

2280 

30000000 

40000000 

50000000 

60000000 

70000000 

80000000 

90000000 

2.6 

2.1 

1.8 

1.5 

1.4 

1.3 

1.2 

2340 

2370 

2410 

2410 

2410 

2410 

2410 

100000000  : 

200000000  : 

300000000 

400000000  i 

500000000  1 

BOOOOOOOO 

700000000 

1,1 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

2440 

2470 

2470 

2470 

2470 

2470 

2470 

800000000 

900000000  1000000000  1100000000  1200000000  1300000000 

1400000000 

0,7 

0.7 

0.7 

0.7 

0.7 

0.7 

0,7 

2470 

2470 

2470 

2470 

2470 

2470 

2470 

1500000000  1600000000  1700000000  1800000000  1900000000  2000000000 

3000000000 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

2470 

2470 

2470 

2470 

2470 

2470 

2470 

4000000000  5000000000  6000000000 

0.7  0,7  0.7 

2470  2470  2470 
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DIM*:(VSIUM  HADI(  S00).TbMP(^0n),BLK(49#2E>),8LKl(49.25)»PTRAN(116). 
1FI  L('116  ) .  rjTPAf\i(  1 1<S  )  ^XMORTI  <  9  n  )  .  TFMP 1  (  9  n)  ,  R  A  D  1  I  (  9  fl  )  ,  ABsMO  ( 1 1 6  )  .  A  PS  I 
2(  11F>  >  »AHS  I  Hll(S),nELTT(9),FTTMs(Q),LEX{22)*OMF(>A(j.  0>*SIGMA(10),TRA 

_ 3N(.116..) _ 

iPARb  =  1 

RFAP  t74,.(PTRAN(  I  ).  I  =  _ _ _ 

PFAD  174,.((JTRAN(  I  ),  1=1 

HFAQ  1  7^.  (  XlvigPT  1  (  I  )  .  1  =  1  .4ft  ) _ 

HF AD  178,  ( TFMPl ( I  >  ,  1=1  , 88  > 

>^FAn  l80,(H/iDri(I),I=1.88) _ 

WFAL)  18?,  (  AHSmu(  I  )  ,  1=1  ,  11(8  > 

WFAP  18?,  ( ANSI  (I  ), 1=1.118) _ 

kFAD  182,  (  AhSI  I  (  I  ) ,  1=1  .11(8  > 

9  PFAp  II  .  M  WIJf'UMCK  1  )  ,  1  =  1,  <  CL  ASSF  n  )  ,  I  =1 , 3)  ) _ 

FFAP  1 8 ft, PFT, cm oh, test, PEI  7,nEl  H , HDl , I M , I N , K K , L 0 N G , L A T 

10  HFAlj  l88.Ylt"LU,ALT.IiI<;T.FNIlF:P.ANSLF,F0CAL,AMUL,AMFL.AMFH,HELHll,CAr 
I^KMFI 

HF  AL'187,  COIN  ,  COK'? _ 

IF  (  cFixiSF  I  TCm  4)14,11 

14  ^^FAP  17b,  ( 1  PX(  I  )  ,  1  =1  ,?0  ) 

DO  13  L  =  1  , 1 

15  08’h0A(L)  =  !l  .  II _ 

11  HFAt  1h9,  (Ulr-LTI  (  I  ),  1=1  ,9) 

HFAP189,  (FT  185(1), 1=1,9) 

S  F  N  9 1  L  1 1>  (1 T  1 

X  =  n  .  P _ 

TsO  .  0 

LFX(?1)  =  0 _ 

LFX(22)  =  0 

K  f  I  =  1 _ 

(>n  TO  413 

12  F'Tp'b  =  FTlMSfKO) 

PFLT  =  Ptl  TI (KO) 

MIJT=F  TIM^:/!1El,  T 

13  T»'A  x  =  (  .  fJ3?*Y  I  FLI'**  .  5) 

IF  (  vlFL''-10  .  )1ft,l7,l7 

18  T^'AX  =  TMAX*1 .3 
17  J  =  1 

im  ?3  K  =  1 ,  ivl  UT 
X  =  X*L1FlT 
y=x/ Tmax 

22  IF(  Y-xmlRTl  (j)  )  24,24. ?3 
■?3  J  =  .J*1 

i>0  TO  22 

“?4  Tt'ip<K)  =  TFMPl(,j) 

?5  H4IJI(K)=HADTI(..i)»yIELT'**.4 
D  W  A  '/F  =  1  .  I)  F  -  ^ 

CnivJ7  =  CUiv*<vMF  I«DIST*MEl  HU 
^COO«  =  AMFl-A>‘F  h 
GO''w.b  =  CUf>/?*U  I  ST^HFL  Hii 
IF(<;FMSt  SWITCH  3)213,27 
?7  IF( ALT)?08,32,?8 
~26  A=(1  . +(UIST/Al  T  )tt*2)**.6 
un  .<U  L  =  1,116 
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_l_AWTnM  Rl  ACKHIIBNI  JR. 


xn  TPAN(L)  =  (  (1.Q-C0N4»(1.0-PTBAN(I  )  )  )  »*A  )  »  (  QTWAN  (  L)  »»Cr)Na  ) 

GO  TO  35 

.^2  S  =  D1ST*AMFL/5.?5 _ 

DO  34  M=l,116 

.■^4  T»Ak(M)s(PTRAN(M  )  »»S  )  » (  QTP  AN  /  M  )  »»CON5) _ 

^5  WFWliMD  B 

Yj  =COSF(  angle  ) _ _ 

Zls. 2245-12 

''1=1.44  5-3 _ 

40  no  47  N=l,ivIDT 

41  BC=.^5.E-6 _ 

DO  46  J!=1,316 

Fl  =?1»TWAM(  Jl  )/(H0»»4tt(?.71B_?ft<H*(X1/  (Bn«TFMP(M)  )  )-!.)> _ 

F|.L(  JI  )  =F  L*  YI/FMUMB#*? 

_ HOshC-t-DwAVfe _ 

46  SUMsSUM  +  FLl-C  JI  )*DWAVE«DELT 

SliMr  =  5Ui‘'i»FiNIIMB«»2»WAUT(M)«»2/(<=i?flf).n»»?)/(AlT«»?-fni‘ST««P) 

47  WRITE  OUTPUT  TAPE  8 » loQ , ( FLL ( K  )  , K  =  1 , ll6 > 

lF(KQ-l)4ll.41Q.4n _ 

410  PRl^T  301,  (PUNiNOC  I  )  ,  1=1.3) 

PP^^T  3  (1 2 _ 

411  PRIM  303,SUf-1.SUMD.FTTME 

2n2  LATA  =  LAT-1 _ 

203  LONGA  =  LONG-l 

213  MT  Qg  (  LQ^G*!  )/2 _ 

214  TS=^ID 

go  N1  s(Hb  r/2.  )/DELZ-t-TS _ 

N?s(KET/2.1'CM0k)/DELZ  +  TS 

N3aTS-(PEr/2.  I/DELZ _ 

C3  =3 ./ (2. /DELZ»*2  +  2./PELR**?fl.71 .4E-3/DELT) 

REWIND  H _ 

NOshOl 

IFISENSP  LIGHT  1  )57.63 _ 

t;7  DO  62  1=1, LONG 

DO  62  J=1,LAT _ 

6l  bL  K(I , J  )=  r 

...6.2  bl:<1  ( I^J)  =  r _ 

6  3  DO  126  t<  =  l,Mf)T 

PFAr.  INHDT  TAPE  H  .  190  .  f  ELI  (  1)  .  I  =  1  , 1  1 6  ) _ 

A  I  =  0  .  C 

AI  I  =  AI _ 

DO  71  J=l,ll6 

a1  =  FLL(J)«DwAVE/RET«(1.-APS^^O(J))»ARsI(J) _ 

A  I =A I +A1 

A2=FLL(J)»UwAVF/CH0R»fl.-ARSM0(J))»(l.-ABSl(J))»AHSII(J) 

71  at  I  =  AI  i+A2 

R  =  F0CAL*RAUI(x-)«AMUL/fDIST#*2  +  ALT*»2)*». 5/5280. 

IF (K-3. 55-4)72,73, 73 

72  A  I sfl ( R/ 3 . 5E-4  1  *»2 
AII=AII»(P/3.5E-4)**2 

73  T  =  T->-DFLT _ ' 

Ll  =  l  .  +R/UELP 

79  I  =  Mr  D _ 

a  1  G  =  ()  .  0 


J=2 
I  J  =  K  I  D 
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R3  IF( I-IJ)89« B4,89 

B4  IF  <  J-1  )H6.85#Q(S _ 

«5  IF( I-(MID-l) )96.96,92 

fl6  IF(  J-JI )ln4,87.B7 _ 

«7  IF( J«H,96»pa 

flQ  IF(  )ino.<?2.lon _ 

R9  IF ( J-1 )90. ino, 90 

Qi.)  IF(I-J-(M1L1-I))l(l8,ln4,lnp _ 

9l  IF< I-IJ)96, ^2.92 

Q2  Ijal _ 

JI=J 

1=1*1 _ 

GO  TO  114 

06  ij=r _ 

JI  =  J 

1  =  1-1 _ 

GO  TO  114 

mo  1.1=1 _ _ 

JI  SJ 

J=J*1 _ 

Go  TO  114 

ln4  Ijal _ 

JI  =J 

J=J-1 _ 

GO  TO  114 

ina  IF(  mj-iMiu-*'!)  )9i,io4,9i _ _ 

114  IF ( J-Ll )115.115,123 

115  IF< I-N2)116.1?1.123 _ 

116  IF(  r-M)117.119,121 

117  IF( T-Ni)l?3,119,119 
119  P1sCAF*aI/1.4F-3 

Go  TO  1124 

mi  PI  sC:mF*a  I  1/3  .4fc-3  ^ 

Go  TO  112<i 

123  PI =0.0  ^ 

1124  lF(SbNSP  SWITCH  5)112';, 124 

1125  PPl'\t  1126,  AI,  All, PI 

1126  FnW^-AT(3H12.4) 

124  IF ( J-1 ) 129, 129, 132 
ip9  \/3  =  --'LK  (  I ,  J>1  ) 

\/c,  =  v/3 
1.^0  SrlCi. 

GO  TO  134 

132  V5  =  PLK  (1,1-1  ) 

V.Ts'-'LK  (  1  ,  J  +  1  )  ■ — ^ — — 

133  Aa=j-1 

- STTOTSTJETR - - - - — 

134  BkO  =  Ci*((wLk(I-i,J)+8!  K(I  +  i,  i)i/D£:L7**?+(V5  +  V3)/DFLR**2-(V5-V3>/2. 
j  /  c/ i’bLH  +  MLK  1  (  !  ,  J)/1.4H-3/!’FLt+p1) 

D  I  U  =  a  f\0-HLK  (  I  ,  J  ) 

If  (■  A  dSF  (  UIlJ  )  -H  TG  )  139,139,1 .3H 
13  9  B I  6  =  A  4  S  F  (  1)  I  D  ) 

'139  c5L-r(  I  ,‘J)  =H  kO 
140  lF(I-2)B3,14l,Ov3 
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1/1?  Dnl43NN=l .  I.aTA _ 

143  HLK(  1,  NN  )  =Hl,K(  ?,  NN  ) 

D014i?NN  =  2.l.ONQA _ 

14  5  BLK<  iMN,LAT  )=BLK(  NN,LAt-1  ) 

D0147N)N  =  l,l.ATA _ 

147  HLK ( LONG, NN ) =BLK ( L0NG_1, NMS 

lF<SbNSfc  TCHl  >?17. 1  49 _ 

jdO  IF <HiG-TEST«BLK( MID, 1) >218,21 8,79 

It;!  DO  153  NN  =  l,LONG _ 

DO  153  NM=1.LAT 

153  HlKI  (  NN,  NIM  )  =bLK  (  NN,  NM  ) _ 

iFISfcNSF  SWITCH  4)150,154 

150  I  =  LEXC 1  ) _ 

J  =  LtX<  2 ) 

L  =  0 _ 

400  L  =  L+1 

SIGMA(L)  =  FXPF(  13  7 . 8712-48650 . 0/309  .  ■‘•8LK  (  I  ,  J  )  1-LOGF  (  DFLT  )  ) 
OMbGAfL)  =  CMEGAIL )+STGMA(L ) 

I  =  LE:X(2«L->1) _ 

152  J  =  LtX(2*L+2> 

IF(I~0)40U,154.40Q _ 

154  IF  (  NO-imOI  )  155. 158, 155 

m5  Nn  =  NiO-t-i _ 

156  Continue 

_ End  f i l  r _ 6 _ 

C  F.  D  F  I  L  c  V 
t  Ni  U  F  I  L  F  H 


sKO  +  1 

IFinbl.TI  (KO)  )?n8.171.12 _ 

153  count  =  0.0 

WCiTb  OUTPUT  TAPE  6  ♦  3  n  4 .  (  (  9UNN0  (  I)  .  I  s  1 .  (  CL  ASSF  (  I  ) ,  I  s  3  ,  3  )  > 

iPAGb  =  Ic'AGE  *  1 

COUNT  =  COUNT  *  1.0 _ 


415  JFK  =  1 
KFJ  =  5 


GO  TO  416 

420  IF(KN  -  l!-n42I  .42?.42? _ 

4?1  KM  =  KN  +  KK 
GO  TO  414 

422  W«ITb  OUTPUT  TAPE  6,3n7 

COUNT  =  COUNT  1.0 _ 

IFICOUNT  -  59 . 0 )422, 423,208 


4?3  WPlTb  OUTPUT  TAPE  6 , 3n9  >  (  Cl.  AgSF  (  I  > .  I  =  1 . 5  ) . . .  - 

NO  =  1 

IFCSENSF  switch  4)434. ts^ _ 

4?4  WRITE  OUTPUT  TAPE  6 , 3 04» ( ( RUMNO ( 1) , I =1 . 3 ) » ( CL ASSF ( 1) . I s 1 , 3 ) ) , I P A SE 

iPAOh  =  IPAOE  1 _ 

WRITE  OUTPUT  TARE  6 , 232»  (  OMEO,  A  (  1)  ,  L  E  X  (  2#  I  - 1 )  ,  LE  X  (  2*  I  )  ,  I  s  1,10) 

WRiTfc  OUTPUT  TAPE  A  .  3nfl.  (  01  ARSF  (  I  )  .  I  =1 . 3  ) _ ' 

GO  TO  156 

495  WRITE  OUTPUT  TAPE  6 . 3fi9.  (  Cl  ASSF  (  I  )  .  1  si .  3  ) _ 

WRITE  OUTPUT  TAPE  6 , 3n4 .  (  ( RUNNO < I ) , I =1 . 3 ) , < CL ASSF ( I ) , I  si , 3 ) ) , I  PAGE 

IPAGE  =  IPAGE  f  1 _ 

COUKTsl. 0 

GO  TO  4lS _ 

4?6  write  OllTPUT  TAPE  6 , 309>  (  Cl  A^SF  (  I  )  ,  1  si ,  3  ) 

WRITE  OUTPUT  TARE  6 , 3  n4 ,  (  (  RUNNi)  (  I  )  .  I  si ,  3  ) .  (  CL  ASSF  (  I  ) ,  I  s  1 , 3  )  ) ,  I  P  AGE 

iPAGE  s  image  4.  1 

CQUNTsl .  0 _ 

GO  TO  417 

1a;5  write  OUTPUT  TAPE  7.1o4,T _ 

DO  168  KNsIm,IN,KK 

WRITE  OUTPUT  TAPE  7,1q6,KN _ 

168  WRITE  Output  tape  7  ,  loo  .  ( RI  K ( KN , NK > , NKsI , L at  ) 

GO  TO  149 _ 

171  IF(SENSE  SWITCH  2)170,9 

170  PAUSE  77777 _ 

172  GO  TO  9 

174  FORPATI l4Fb  .  3  ) _ 

175  FORmAT(20I3) 

176  FORMAT! 7F10 .8  ) _ 

178  FORMAT! 7F1 0 .1 ) 

180  F0RMAT(16F4.1) _ 

1P2  FORMAT! 7F 10 ,6  ) 

^  F0RmAT!5F6.6.6I4) _ 

It  F0RmAT!2F4,3) 

Iflb  F0RmAT!F  9.2.F8.3,F9.0,F5.2,F6.3,F7.3»F6.3#5F4.3) 

1Q9  FORMAT! 7F10 . 7  ) 

190  FORMAT! 3E14. 5) 

191  FORMAT! 5tl4. 7) 

192  FoRmAT!30M,  END  OF  RELAXATION  FOR  TIME  s  F7.6) 

194  F0RMAT!1H  8ri  TIME  =  F9.7,8H  R  =  E12.5) 

196  F0RMAT!6h  I  =  12) 

198  FORMAT! 5n  5F15.61 

200  FORMAT!////) 

208  PAUSE  123456  ~ 

GO  TO  27 

217  IFlf.U-N01)l49,16$.l49 
210  IF(SENSE  switch  1)219,151 


219  IF  !  NO-NO jT)  151, 220, 151 

220  write  output  tape  7»l92,T 
Go  to  151 

413  WRITE  OUTPUT  TAPE  6 , 3 n4 , ! ( RUNNO ( I ) , I s 1 , 3  ) , < CL ASSF !  I ) , I s 1 , 3 ) ) , I P AGF 
IPAGE  s  IPAGE  1  ^  ^ 

_  WRITt  OUTPUT  TAPE  6,2n0 

224  WRITE  OUTPUT  TAPE  6 . 228 , Y I FL 0 , ALT , D I ST , F NUMB , ANGLE , FOC AL . AMUL , AMFL 

1,  AMFH,RELHu,WMFI,CAF,PET,CHOR,TEsT,nELZ,DELR,N01, IM, I N, KK, LONG, LAT 

2, DEl.TI!i),DELTI!2)/l)El  TI13)»DEI.TI!4).I)ELT1(5),DELTI(6),DELT117),DE 
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write  output  tape  6.200 

WRITE  0  U  T  P  UT  TAPP  6 . 2.-^  0  .  (  F  T  I  MS  M  )  .  I  =1 . 9  ) 


WRITE  OUTPUT  TAPfc  6 , 3  0  (  Cl,.  Ac;SP  (I  )  ,  I  si ,  3  ) 

SLlMsO  . 


an  TO  12 

MAT(  RM'YlPiri  =  FR.2/9H  AlT  =  F10.3/9H  DIST  =  FIO.O/QH  FWII 


=  Fin. 2/RH  angle  =  F10.3/9H  FOCAL  =  F10.5/9H  AMIII.  s  F10.3/9H  AMF 

I  Hll  =  Fln.3/<}H  WMFT  =  FlO-.X/QH  C 


3AF  =  F10.3/QH  RET  s  F10.6/OH  ChOR  =  F10.6/9H  TEST  =  F10.6/9H 
d  nPl  7  =  Flri.(«,/9H  DFID  =  Mf)l  =  IIO/QH  IM  s  Tltl/OH  t 


5jN  =  I10/9H  KK  =  I10/9H  LONG  =  I10/9H  LAT  =  I10/12H  DELTK 

=  rin.7/1?H  DFlTIfX)  =  Fln.7/15H  DPI 


=  F10.7/12H  DELT1(5)  =  <^10.7/1?H  DFLTI(6)  s  F10.7/1?H  DELTII7)  = 
ln.7/12H  DElTIlP)  =  Fin. 7) 


230  F0Hf^AT(l2H  FTIMSd)  =  F10.7/12H  FTIMS(2)  =  F10.7/12H  FTIMS<3)  : 
1F10.7/12H  FTIMS(4)  =  Fli).7/1?H  FTIMS(5)  =  FI 


?F10.7/12h  FTIMS(7)  =  FlO. 7/120  FTlMS(fl)  s  FlQ.?) 

FnRMAT(23ri  nMEGA  T  J/(E12,5.3H  1 3 . 3H  13/)) _ 


300  format ( 3A6, 3A6  ) 
nl  FORMATdHI  40X,3A6) 


302  F0RMAT( 28X, 3hSUM. 26X,4HSUMn,25X, 5HFTIM6 ) 
0  3  (■0RMAT(24X.El?,b,t8X,Pl2.5,lflX.El2.5) 


304  f-ORMATClril  1  0  X  ,  3  A6 , 1 0  X  ,  3  A6 , 21  X .  6H  PAGE  14) 
5  FORMAT!///////////  40X.3AA) 


307  FORMAT!/) 


309  format ! 40X, 3A6 ) 

ENU  <0.1.  (1.0. 1,0. (1,1, 0.0. 0,0. 0,0,0) 


END  (0,1, 0,1,1) 


♦N  END  TAPE 


TABLE  3 


PTRAN  (atmospheric  transmission  less  water  vapor)  (from  350  to  1500  microns  in 
increments  of  10  microns)  (  x  lO"*) . 


551 

589 

623 

656 

685 

712 

746 

758 

776 

795 

812 

827 

839 

852 

863 

874 

882 

891 

898 

906 

913 

918 

923 

929 

934 

938 

941 

945 

949 

952 

955 

957 

959 

962 

964 

966 

968 

969 

970 

972 

974 

976 

977 

978 

979 

980 

981 

981 

982 

983 

984 

984 

985 

986 

987 

988 

988 

989 

989 

990 

990 

990 

991 

991 

992 

992 

992 

993 

993 

993 

994 

994 

994 

994 

995 

995 

995 

995 

995 

995 

995 

996 

996 

996 

996 

996 

996 

996 

996 

996 

997 

997 

997 

997 

997 

997 

997 

997 

997 

997 

998 

998 

998 

998 

998 

998 

998 

998 

998 

998 

998 

998 

998 

998 

998 

998 

QTRAN 

(water  vapor 

transmissivity) 

(from 

350  to 

1500 

microns  in 

1,?^  “■  ; 

■■■.i.ts  of. 

10  microns) 

(  X  10  ). 

926 

934 

940 

945 

949 

951 

953 

955 

957 

959 

981 

953 

964 

368 

967 

968 

970 

971 

972 

971 

971 

970 

970 

971 

9"^’ 2 

973 

974 

975 

976 

977 

978 

979 

980 

981 

981 

982 

982 

983 

983 

983 

984 

984 

985 

985 

985 

905 

985 

986 

986 

9S6 

986 

986 

986 

988 

986 

986 

987 

987 

987 

987 

987 

987 

98? 

987 

987 

387 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

387 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

SS7 

987 

9?  7 

337 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

987 

980 

987 

987 

987 

987 

987 

987 

3 

Absorption  coefficient  in  ocular  media  x  10  .  (ABSMO) 


1000000 

1000000 

1000000 

1000000 

1000000 

1000000 

0900000 

0620000 

0575000 

0500000 

0425000 

0333000 

0285000 

0240000 

0223000 

C 185000 

0160000 

0130000 

0115000 

0100000 

0090000 

0075000 

0060000 

0050000 

0037000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033OC0 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0033000 

0035000 

0040000 

0050000 

0060000 

0075000 

0088000 

0110000 

0150000 

0290000 

0480000 

0534000 

0575000 

0582000 

0582000 

0582000 

0550000 

0555000 

0475000 

0400000 

0300000 

0250000 

0210000 

0200000 

0210000 

0225000 

0260000 

0290000 

0360000 

0410000 

0680000 

0810000 

0880000 

0900000 

0920000 

0930000 

0940000 

0950000 

0940000 

0930000 

0920000 

0910000 

0905000 

0900000 

0905000 

0910000 

0920000 

0940000 

0970000 

0980000 

0985000 

0990000 

0993000 

0995000 

1000000 

1000000 

1000000 

1000000 

1000000 

1000000 

1000000 

1000000 

1000000 

1000000 
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TABLE  4 


Reduced  height,  hi  and  decimal  percent,  of  th'»  atmosphere  for  a  station 
that  is  k  km  above  sea  level.  (First  three  columns  of  this  table  from 
Kulper,  G.P. ,  The  Atmosphere  of  the  Earth  and  Planets,  1952.) 


k 

^1 

(Winter) 

h.  (Summer) 

Decimal  Percent  Decimal  Pen 

(km) 

(km) 

(km) 

(Winter) 

(Summer) 

0 

8.03 

8.05 

1.00 

1.00 

1 

7.09 

7.15 

0.884 

0.890 

2 

6.26 

6.34 

0.780 

0.787 

3 

5.51 

5.61 

0.686 

0.897 

4 

4.83 

4.95 

0.602 

0.615 

5 

4.23 

4.36 

0.526 

0.541 

6 

3.68 

3.82 

0.45S 

0.475 

7 

3.20 

3.34 

0.399 

0.415 

8 

2.76 

2.91 

0.344 

0.361 

9 

2.38 

2.53 

0.296 

0.314 

10 

2.04 

2.18 

0.254 

0.281 

12 

1.45 

1.61 

0.180 

0.200 

14 

1.09 

1.19 

0.136 

0.128 

16 

0.79 

0.88 

0.098 

0.011 

18 

0.58 

0.65 

0.072 

0.081 

20 

0.42 

0.48 

0.053 

0.060 

25 

0.19 

0.22 

0.024 

0.027 

30 

0.09 

0.11 

0.011 

0.014 

35 

0.04 

0.05 

0.005 

0.006 

40 

0.02 

0.02 

0.003 

0.003 

The  above 

table  is 

used  to  find  constants  and  Mjj 

(designated  as 

AMFL  and  AMFH  in  program) 
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TABLE  5 
(Read  Across) 


Absorption  coefficient  in  PE  x  10^.  (ABSI) 


0853000 

0825000 

0805000 

0775000 

0750000 

0730000 

0723000 

0717000 

0714000 

0712000 

0710000 

0700000 

0687000 

0675000 

0665000 

0650000 

0625000 

0600000 

0575000 

0550000 

0530000 

0505000 

0480000 

0457000 

0435000 

0420000 

0400000 

0385000 

0370000 

0355000 

0343000 

0333000 

0315000 

0300000 

0287000 

0273000 

0260000 

0250000 

0240000 

0225000 

0215000 

0205000 

0195000 

0187000 

0180000 

0175000 

0168000 

0162000 

0157000 

0152000 

0150000 

0145000 

0140000 

0137000 

0130000 

0125000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0120000 

0110000 

0110000 

OilOOOO 

0110000 

0110000 

012000C 

0125000 

0134000 

0140000 

0145000 

0140000 

0137000 

0135000 

0127000 

0120000 

0110000 

0107000 

0102000 

0100000 

0100000 

0113000 

0125000 

0140000 

0150000 

0160000 

0175000 

0190000 

0200000 

0215000 

0230000 

0350000 

0510000 

0675000 

0840000 

0975000 

Absorption  coefficient  in 

choroid  x 

10^.  (ABSII) 

0923000 

0923000 

0923000 

0920000 

0915000 

0910000 

0905000 

0900000 

0890000 

0882000 

0875000 

0860000 

0890000 

0820000 

0795000 

0775000 

0765000 

0760000 

0750000 

0745000 

0740000 

0733000 

0720000 

0705000 

0685000 

0660000 

0642000 

0625000 

0605000 

0583000 

0555000 

0547000 

0540000 

0533000 

0527000 

0520000 

0500000 

0475000 

0450000 

0417000 

0390000 

0375000 

0370000 

0360000 

0350000 

0345000 

0340000 

0330000 

0323000 

0317000 

0310000 

0295000 

0283000 

0270000 

0255000 

0245000 

0240000 

0230000 

0227000 

0224000 

0222000 

0230000 

0240000 

0250000 

0257000 

0270000 

0265000 

0260000 

0255000 

0245000 

0230000 

0220000 

0210000 

0197000 

0187000 

0180000 

0175000 

0170000 

0170000 

0170000 

0175000 

0195000 

0225000 

0253000 

0285000 

0315000 

0325000 

0334000 

0340000 

0340000 

0340000 

0340000 

0335000 

0330000 

0320000 

0315000 

0300000 

0295000 

0285000 

0280000 

0275000 

0303000 

0340000 

0390000 

0450000 

0520000 

0690000 

0795000 

0870000 

0935000 

0975000 

0975000 

0975000 

0975000 

0975000 

0975000 
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TABLE  6 


RUN  G  4 


UNCLASSIFIED 


YIELD  » 

2.50 

ALT 

0. 

DIST  « 

1. 

FNUMB  = 

8.50 

ANGLE  = 

0. 

FOCAL  = 

1.70000 

AMUL  = 

1.000 

AMFL  = 

1.000 

AMFH  * 

0. 

RELHU  = 

0.500 

WMF  I.  = 

0. 

CAF 

1.000 

RET 

0.001 000 

CHOR  = 

0.010000 

TEST  = 

0.010000 

DELZ  * 

0.001000 

DELR  » 

0.001000 

NOl 

1 

IM 

1 

IN 

49 

KK 

3 

LONG  > 

49 

LAT  a 

25 

DELTK  I ) 

«  0.0000020 

DELTI(2) 

>  0.0000200 

DELTIO) 

s  0.0002000 

DELT1(4) 

>  0.0020000 

DELTKS) 

>  0.0200000 

CELT  I (6) 

«  0.2000000 

DELTII7I 

»  -0. 

DELTI(0» 

»  -0. 

FTIMSI  1) 

X 

0.0000100 

FTIMS(2) 

M 

0.0001000 

FTIMS(3) 

M 

0.0010000 

FTIMS14) 

X 

0.0100000 

FTIMSC5) 

X 

0.1000000 

FTIMSC6) 

X 

0.5000000 

FTIMSI 7) 

X 

-0. 

FTIMS(8) 

=: 

-0. 

35 


RUN  G  4 

0.0000020  R 


UNCLASSIFIED 

0.78966E-02 


TIME  = 

I  =  1 

0. 

0. 

0. 

0. 

0. 

I  =  4 

0. 

0. 

0. 

0. 

0. 

I  =  7 

0. 

0. 

0. 

0. 

0. 

I  =  10 
0. 15304E-35 
0.  14246E-35 

0. 

0. 

0. 

I  =  13 
0.70935E-28 
0. 67053 E-20 
0.78622E-34 

0. 

0. 

I  *  16 
0.32881E-20 
0.31567E-20 
0. 13466E-26 
0. 13594E-37 
0. 

I  =  19 
0. 15237E-12 
0. 14864E-12 
0.42292E-19 
0. 14341E-30 
0. 

I  *  22 
0.70602E-05 
0.70003E-05 
0.87201E-12 
0.58383E-24 
0. 18191E-36 
I  »  25 
0.91084E  00 
0.91085E  00 
0. 16363E-07 
0.’22167E»20 
0.31787E-33 
I  =  28 
0.67765E-01 
0.57765E-01 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0.  14907E-35 
0.82105E-37 

0. 

0. 

0. 

0.69320E-28 

0.27577E-29 

0.69278E-36 

0. 

0. 

0.32369E-20 
0.77623E-22 
0.  17209E-28 

0. 

0. 

0. 15055E-12 
0.1 1958E-14 
0.35648E-21 
0.45100E-33 

0. 

0.7031 3E-05 
0.56322E-07 
0. 39553E-14 
0. 17899E-26 

0. 

0. 91084E  00 
0.91085E  00 
0.43379E-10 
0.59640E-23 
0.06180E'36 

0.57765E-0I 

0.57765E-01 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0.  152  05E-35 
0.  14196E-35 

0. 

0. 

0. 

0.70493E-28 

0.66837E-28 

0-22892E-38 

0. 

0. 

0.32675E-20 
0.31471E-20 
0.598 25E-31 

0. 

0. 

0. 15143E-12 
0. 14821E-12 
0. 12240E-23 
0.22471E-35 

0. 

0.70163E-05 
0.69803E-05 
0.  131 19E-16 
0.61405E-29 

0. 

0.9  10B5E  00 
0.  90817E  00 
0. 1 1545E-12 
0. 16076E-25 
0.23385E-38 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0.  14445E-35 
0.33098E-37 

0. 

0. 

0. 

0.67830E-28 
0.  10394E-29 

0. 

0. 

0. 

0.31785E-20 

0.24463E-22 

0.54929E-33 

0. 

0. 

0. 14931E-12 
0.76846E-15 
0.79536E-26 
0.69132E-38 

0. 

0.69990E-05 
0.36193E-07 
0.51287E-19 
0. 18477E-31 

0. 

0.91085E  00 
0.23645E-02 
0.30831E-15 
0.43409E-28 

0. 


0.57765E-01  0.57764E-01 

0.57595E-01  0.15034E-03 

36  UNCLASSIFIED 


0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0.97883E-37 

0. 

0. 

0. 

0. 

0.34972E-29 
0.  12683E-31 

0. 

0. 

0. 

0.11227E-21 
0.34140E-24 
0.  18024E-35 

0. 

0. 

0.28225E-14 
0. 11057E-16 
0.25851E-28 

0. 

0. 

0.55506E-07 
0.23683E--09 
0. 16185E-2I 
0.61248E-34 

0. 

0.91085E  00 
0.62020E-05 
0.82567E-18 
0.11739E-30 
0. 

0.57765E-01 

0.39536E-06 


RUN  G  4 


UNCLASSIFIED 
0. 74 163E-14 
0. 10457E-26 


0. 10458E-08 
0. 14347E-21 
0.2083JE-34 
I  =  31 
0.57766E-0  1 
0.57766E-01 
0. 10457E-08 
0. 14345E-21 
0.20822E-34 
I  =  34 
0.57765E-0 1 
0.57765E-01 
0.  10343E-08 
0.  13887E-2 1 
0.  19797E-34 
I  =  37 
0.44661E-06 
0.4466 lE-06 
0.23747E-13 
0.74  764E-26 
0.23092E-38 
I  =  40 
0.96400E-14 
0.96402E-14 
0.46929E-20 
0.28261E-32 
0. 

I  =  43 
0.20808E-2 1 
0.20807E-21 
0.24915E-27 

0. 

0. 

I  =  46 
0.44905E-29 
0. 44907 E-29 
0. 13769E-34 

0. 

0. 

I  =  49 
0.34900E-34 
0.34902E-34 

0. 

0. 

0. 


0.27794E-n 
0.38698E-24 
0.  5661  6E-37 

0.57766E-01 
0.57765E-0  1 
0.27793E-1  1 
0.38690E-24 
0.565a8E-37 

0.57765E-01 
0.57765E-01 
0.27340E-1 1 
0.373608-24 
0. 53524E-37 

0.4466  lE-06 
0.  44659E-06 
0.6288  IE-  16 
0.25915E-28 

0. 

0.96402E-  14 
0.96389E-14 
0.  19944E-22 
0.  1  1873E-34 

0. 

0.20808E-21 

0.208018-21 

0.127018-29 

0. 

0. 

0.449078-29 
0.448858-29 
0. 923338-37 

0. 

0. 

0.349028-34 

0.348808-34 

0. 

0. 

0. 


0. 

0.57766E-0 1 
0. 57595E-0 1 
0. 741578-14 
0. 10455E-26 

0. 

0. 57765E-01 
0.57595E-01 
0.72756E-14 
0. 10042E-26 

0. 

0.4466 lE-06 
0.442688-06 
0. 16748E-18 
0. 77469E-31 

0. 

0.96402E-14 

0.94729E-14 

0.776708-25 

0.349978-37 

0. 

0.208088-21 

0.202728-21 

0.576888-32 

0. 

0. 

0.449078-29 

0.433898-29 

0. 

0. 

0. 

0.349028-34 

0.33539E-34 

0. 

0. 

0. 


0.  198558-16 
0.28307E-29 

0. 

0.577668-0 1 
0. 150348-03 
0. 198538-16 
0.282988-29 

0. 


0.577658-01 
0. 149938-03 
0.193738-16 
0.271138-29 

0. 

0.446618-06 

0.343288-08 

0.594828-21 

0.25175E-33 

0. 

0.964028-14 
0.  146278-15 
0.279018-27 

0. 

0. 

0.208088-21 

0.466798-23 

0.240788-34 

0. 

0. 

0.449078-29 

0.132638-30 

0. 

0. 

0. 

0.349028-34 

0.119078-35 

0. 

0. 

0. 


0.533078-19 

0.767388-32 

0. 

0.577668-01 

0.395358-06 

0.533008-19 

0.767108-32 

0. 

0.577658-01 
0.392108-06 
0.51876E-19 
0.731  148-32 
0. 

0.446618-86 
0.90023E-1 1 
0.239228-23 
0.73628E-36 

0. 


0.964028-14 

0.955528-18 

0.988578-30 

0. 

0. 


0.208078-21 

0.406098-25 

0.948768-37 

0. 

0. 

0.449078-29 
0.  164148-32 

0. 

0. 

0. 

0.349028-34 
0.  17042E-37 

0. 

0. 

0. 
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UNCLASSIFIED 


RUN  G  4 

0.0000300  R  =  0.24619E-01 

0. 14631E-33 


TIME  = 

I  =  1 

0.15589E-33 
0. 96530E-34 
0.'37374E-34 
0.442nE-34 
0.  10443E-34 
I  =  4 

0.21866E-30 
0.  I3821E-30 
0.49088E-31 
0.65321E-31 
0. 99894 E-32 
I  =  7 

0. 1 18 14E-25 
0.77880E-26 
0.23219E-26 
0.38957E-26 
0.43458E-27 
I  =  10 
0.64266E-21 
0.44185E-21 
0. 10120E-21 
0.23696E-21 
0.27018E-22 
I  =  13 
0.34637E-16 
0.25468E-16 
0.37083E-17 
0. 14616E-16 
0.16796E-17 
I  =  16 
0.18752E-11 
0. 14856E-1 1 
0. 12517E-12 
0.90133E-12 
0. 10442E-12 
I  =  19 
0. 10081E-06 
0.88533E-07 
0. 85464 E-0a 
0.55575E-07 
0.64909E-08 
I  =  22 
0.53874E-02 
0-52912E-02 
0.41343E-02 
0.34262E-02 
0.33329E-02 
I  =  25 
0.67417E  01 
0.67422E  01 
0.58470E  01 
0.S0096E  01 
0.500956  01 
I  =  28 
0.44162E-00 
0.44156E-00 


0. 15540E-33 
0. 66356E-34 
0. 54505E-34 
0. 18429E-34 
0.41069E-34 

0.2186  lE-30 
0.89774E-31 
0.79440E-31 
0.22044E-31 
0.62222E-31 

0.  1  186  3E-25 
0.45262E-26 
0.46269E-26 
0.80515E-27 
0.38396E-26 

0.64396E-21 

0.21990E-21 

0.27406E-21 

0.26981E-22 

0.23687E-21 

0.  34781E-  16 
0. 10089E-16 
0. 16552E-16 
0. 16773E-17 
0.  14612E-16 

0.18777E-11 
0.41329E-12 
0.  101826-n 

0.  10427E-12 
0.901  18E-12 

0.  10095E-06 
0. 12793E-07 
0.62773E-07 
0.648 17E-08 
0.55574E-07 

0.53849E-02 

0.51842E-02 

0.38386E-02 

0.33323E-02 

0.34266E-02 

0.67419E  01 
0.674  12E  01 
0.54453E  01 
0.50095E  01 
0.50095E  01 

0.44163E-00 
0.  441  56E-00 


0-82898E-34 

0.29132E-34 

0.42713E-34 

0.75889E-35 

0.20672E-30 
0. 1 1979E-30 
0.37504E-31 
0-63635E-31 
0. 76589E-32 

0.  1  1324E-25 
0.68653E-26 
0. 16906E-26 
0.38470E-26 
0.47020E-27 

0.62061E-21 
0.39835E-21 
0.6679 lE-22 
0.23680E-21 
0.28856E-22 

0.33885E-16 
0.23496E-16 
0. 18295E-17 
0. 14608E-16 
0.  17703E-1  7 

0.  1848  lE-1 1 
0. 14098E-1 1 
0. 11366E-12 
0.90098E-12 
0. 10856E-12 

0. 10028E-06 
0.86317E-07 
0. 72698E-08 
0.55564E-07 
0. 66546E-08 

0.53773E-02 
0.52160E-02 
0. 36953E-02 
0-34261E-02 
0.3337 lE-02 

0.67422E  01 
0.67065E  01 
0.54222E  01 
0.50095E  01 
0. 50095E  01 


0. 12610E-33 
0.50409E-34 
0.48003E-34 
0.  14397E-34 
0.645 15E-34 

0. 17868E-30 
0.67231E-31 
0.70412E-31 
0.  15984E-31 
0.94740E-31 

0.98147E-26 

0.32907E-26 

0.41476E-26 

0.43440E-27 

0.54968E-26 

0.54231E~21 
0.  I5241E-21 
0.24885E-21 
0.26997E-22 
0.31830E-21 

0.29958E-16 
0.63865E-17 
0. 15240E-16 
0. 16784E-17 
0. 18396E-16 

0.  16721E-1  1 
0. 21038 E-12 
0.93977E-12 
0. 10434E-12 
0. 10612E-1  1 

0.935826-07 
0. 10054E-07 
0.57940E-07 
0.64864E-08 
0.61 1 1 lE-07 

0.52854E-02 

0.42339E-02 

0.35605E-02 

0.33326E-02 

0.35133E-02 

0.67422E  01 
0.58941E  01 
0.51675E  01 
0.50095E  01 
0.500956  01 


0. 96962 E-34 
0.64494E-'34 
0.23054E-34 
0.41690E-34 

0. 

0. 13489E-30 
0.933 72 E-31 
0.288 29 E-31 
0.62641E-31 
0.94740E~31 

0-71518E-26 
0.53751E-26 
0.  12043E-26 
0.38390E-26 
0.54968E-26 

0.37432E-21 

0.31416E-21 

0.39355E-22 

0.23684E-21 

0.31830E-21 

0.  19281E-16 
0. 18703E-16 
0.  17178E-17 
0. 14610E-16 
0.  18396E-16 

0.96772E-12 
0. 1 1360E-1 1 
0. 10677E-12 
0.90109E-12 
0. 10612E-11 

0.46709E-07 
0. 70003E-07 
0.66369E-08 
0.55569E-07 
0.6111 lE-07 

0.51831E-02 
0.42685E-02 
0.  334nE-02 
0.34264E-02 
0.35133E-02 

0.67423E  01 
0.58633E  01 
0.50141E  01 
0.50095E  01 
0.50095E  01 

0.44163E-00 

0.38380E-00 


0.44162E-00  0.44153E-00 

0.439 18E-00  0.38593E-00 

38  UNCLASSIFIED 


RUN  G  4 
0.38277E-00 
0. 32687E-00 
0.32693E-00 


UNCLASSIFIED 


I  =  31 
0.44153E-00 
0.  441  52  E-00 
0.’38267E-00 
0.32686E-00 
0. 32685E-00 
I  »  34 
0.44 1 18E-00 
0.44 1  19E-00 
0.38 2 40 E-00 
0.32664  E-00 
0.32664E-00 
I  »  37 
0. 34370E-03 
0.34368E-03 
0.27599E-03 
0.21205E-03 
0.2181 5E-03 
I  =  40 
0.64635E-08 
0. 64638 E-08 
0.51201E-08 
0.41346E-09 
0.35380E-08 
I  =  43 
0. 12062E-12 
0. 12055E-12 
0.94815e-13 
0. 15071E-13 
0.57372E-13 
I  =  46 
0.22453E-17 
0. 22450 E-1 7 
0- 17527E-17 
0.38141E-ia 
0.93235E-1 8 
I  =  49 
0. 16154 E-20 
0. 16153E-20 
0. 12526E-20 
0.31 1 15E-21 
0.61836E-21 


0. 35620E-00 
0.  32693E-00 
0.32686E-00 

0.44153E-00 
0.4414  IE-00 
0.356iaE-00 
0.32685E-00 
0. 32605E-00 

0.4411 9E-00 
0.44  1  12E-00 
0. 35594E-00 
0. 32664E-00 
0.32663E-00 

0. 34370E-03 
0. 34342E-03 
0.25018E-03 
0.218126-03 
0.21208E-03 

0.64645E-08 
0.6449  3E-0a 
0.  46634E-08 
0.35373E-0a 
0.41337E-09 

0.  12062E-12 
0.  12005E-12 
0.86345E-13 
0.577616-13 
0.66497E-14 

0.22463E-17 
0.22310E-17 
0. 15913E-17 
0.962a2E-18 
0.107006-18 

0.161656-20 
0.160306-20 
0. 1 1343E-20 
0.649726-2  1 
0. 128866-21 


0.  3547 lE-00 
0.  32686E-00 
0.32693E-00 

0.44153E-00 
0. 43913E-00 
0.35462E-00 
0.32685E-00 
0.32685E-00 

0.44  1  19E-00 
0.43881E-00 
0. 35439E-00 
0.  32663E-00 
0.32664E-00 

0.34370E-03 
0. 33873E-03 
0.24196E-03 
0.2  1205E-03 
0.21791E-03 

0. 64648E-08 
0.62976E-08 
0.42322E-08 
0.41280E-09 
0.352816-08 

0. 12061E-12 
0. 1 1635E-12 
0.74934E-13 
0.94395E-14 
0. 571 17E-13 

0.22467E-17 
0.21507E-17 
0- 13393E-17 
0-28819E-18 
0.92459E-18 

0. 16169E-20 
0. 1 5409E-20 
0.93961E-21 
0.248586-2  1 
0.60701E-21 


0.33773E-00 

0.32693E-00 

0.32686E-00 

0.44 152E-00 
0.38582E-00 
0.33771E-00 
0.32685E-00 
0.32685E-00 

0.44 120E-00 
0.38555E-00 
0.33749E-00 
0. 32664E-00 
0.32663E-00 

0.34369E-03 
0.28280E-03 
0.22093E-03 
0.218  14E-03 
0.20208E-03 

0.64648E-08 
0.53552E-08 
0.94136E-09 
0.353  77E-08 
0.25265E-09 

0. 12060E-12 
0. 10079E- 12 
0.24442E-13 
0.57361E-13 
0.40674E-14 

0.22467E-17 
0. 18902E-17 
0.55349E-18 
0.94315E-18 
0.65479E-19 

0. 161 70E-20 
0. 13652E-20 
0.43388E-21 
0.631  12E-21 
0.42790E-22 


0.327256-00 
0. 326866-00 
0.326866-00 

0.44 1526-00 
0.383776-00 
0.327176-00 
0.326856-00 
0.326856-00 

0.441206-00 

0.383496-00 

0.326966-00 

0.326636-00 

0.326636-00 

0.343696-03 

0.278446-03 

0.21871E-03 

0.212076-03 

0.202086-03 

0.646476-08 
0.520276-08 
0.356436-08 
0.4131 16-09 
0.252656-09 

0.  120596-12 
0.969776-13 
0.591616-13 
0.664536-14 
0.406746-14 

0.224656-17 
0.  180216-17 
0. 100296-17 
0.197376-18 
0.65479E-19 

0. 16168E-20 
0.  129516-20 
0.683646-21 
0. 18837E-21 

0. 


39 


unclassified 


KUN  O  4 

TIME  =  0.0008900  R  *  0. 

I  »  1 

0.25187E-12  0.24429E-12 

0.  18407E--12  0.16183E-12 

0.81461E-13  0.67398E-13 

0.32252E~13  0.26958E-13 

0.14250?£-13  0.  15255E-I3 

I  =  4 

0.55691E-n  0.54417E-11 

0.4l674E-n  0.36724E-n 

0.  18607E-n  0.1  5434E-11 

0.74970E-12  0.62173E-12 

0.30254E-12  0.40705E-12 

I  *  7 

0.651 16E-09  0.64692E-09 

0.51094E-09  0.45217E-09 

0.23214E--09  0.  19386E-09 

0.98730E-10  0.  76391E-10 

0.47880E-10  0.80763E-10 

I  =  10 

0.72859E-07  0.71679E-07 

0.  57006E-07  0.50648E-07 

0.26444E-07  0.22530E-07 

0.13706E-07  0.926 12E-08 

0.93802E-0a  0.14408E-07 

I  =  13 

0.67a20E-05  0.67603E-05 

0.55932E-05  0.50169E-05 

0.26910E-05  0.25033E-05 

0.21543E-05  0.16123E-05 

0.16223E-05  0.221 16E-05 

I  =  16 

0.53999E-03  0.53643E-03 

0.45654E-03  0.41579E-03 

0.25608E-03  0.288656-03 

0.26307E-03  0.2241  lE-03 

0.22428E-03  0.26491E-03 

I  *  19 

0.30352E-01  0.30348E-01 

0.27509E-01  0.25498E-01 

0.22244E-01  0.22419E-01 

0.20451E-01  0.19513E-01 

0.19470E-01  0.20347E-01 

I  =  22 

0.93170E  00  0.93171E  00 

0.91924E  00  0.90679E  00 

0.79068E  00  0.76472E  00 

0.69955E  00  0.69674E  00 

0.69498E  00  0.69426E  00 

I  =  25 

0.41601E  01  0.41624E  01 

0.41305E  01  0.40767E  01 

0.35977E  01  0.34942E  01 

0.32285E  01  0.32158E  01 

0.32105E  01  0.32116E  01 

I  *  28 

0.79875E  00  0.79750E  00 

0.78870E  00  0.77800E  00 


8037E-01 

0.23456E-12  0.22136E-12 

0.13920E-12  0.11765E-12 

0.55734E-13  0.46262E-13 

0.23051E-13  0.19221E-13 

0. 14220E-13  0.26076E-13 

0.52548E-11  0.49797E-11 

0.31645E-11  0.26785E-11 

0.12782E-11  0.10655E-11 

0.54608E-12  0.43819E-12 

0.31921E-12  0.77718E-12 

0.63 185  E-09  0.60  346E-09 

0.39102E-09  0.33193E-09 

0. 16056E-09  0.13586E-09 

0.78146E-10  0.47046E-10 

0.63883E-10  0.13771E-09 

0.69790E-07  0.66797E-07 

0.43974E-07  0.37463E-07 

0.  18342E-07  0.  16582E-07 

0.13516E-07  0.91956E-08 

0-11648E-07  0.21785E-07 

0.66527E-05  0.64309E-05 

0.44026E-05  0.37771E-05 

0.17532E-05  0.22375E-05 

0.21260E-05  0.16027E-05 

0.18663E-05  0.29360E-05 

0.52722E-03  0.51328E-03 

0.37418E-03  0.31915E-03 

0.24061E-03  0.27262E-03 

0.26015E-03  0.22308E-03 

0.24000E-03  0.30855E-03 

0.30125E-01  0.29695E-01 

0.25706E-01  0.23799E-01 

0.20892E-01  0.21160E-01 

0.20269E-01  0.19448E-01 

0.19828E-01  0.21319E-01 

0.92995E  00  0.92486E  00 

0.  8789I*E  00  0.84459E  00 

0.74318E  00  0.72236E  00 

0.69491E  00  0.69514E  00 

0.69477E  00  0.69346E  00 

0.41585E  01  0.41582E  01 

0.39702E  01  0.38098E  01 

0.34041E  01  0.33210E  01 

0.32129E  01  0.32108E  01 

0.321 19E  01  0.32136E  01 

0.79518E  00  0.79349E  00 

0.75658E  00  0.72724E  00 

UNCLASSIFIED 


0.20434E-12 
0.  98242 E- 13 
0.38438E-13 
0. 17191E-13 

0. 

0.46127E-] 1 
0.22405E-1.1 
0.88499E-12 
0.42438E-12 
0.77718E-12 

0.56256E-09 
0.27872E-09 
0.111326-09 
0.74149E-10 
0.  13771E-0,9 

0.625nE-07 
0. 31689E-07 
0. 12221E-07 
0. 13538E-07 
0.21785E-07 

0.60770E-05 
0.32725E-05 
0.  16536E-05 
0.21273E-05 
0.29360E-05 

0.49040E-03 
0. 30890E-03 
0.22919E-03 
0.26009E-03 
0.30855E-03 

0.289446-01 
0.23915E-01 
0.  19895E-01 
0.20254E-01 
0.21319E-01 

0.92591E  00 
0.81447E  00 
0.70815E  00 
0.69442E  00 
0.69346E  00 

0.41495E  01 
0.36942E  01 
0.32582E  01 
0.32115E  01 
0.32136E  01 

0.79289E  00 
0.70320E  00 


40 


RUN  G  4 
0.68347E  00 
0.60869E  00 
0.60489E  00 


UNCLASSIFIED 
0.64465E  00 
0.60494E  00 
0.60433E  00 


I  =  31 
0.50161E  00 
0.4959  lE-00 
0.43165E-00 
0. 38586E-00 
0.38391E-00 
I  =  34 
0.42891E-00 
0.42583E-00 
0.37 130E-00 
0.33295E-00 
0.331 14E-00 
I  =  37 
0.66464E-0  1 
0.65720E-01 
0.56573E-01 
0.  49802E-0  1 
0.49515E-01 
I  =  40 
0.20755E-02 
0.20591E-02 
0. 17631E-02 
0. 13271E-02 
0. 13655E-02 
I  =  43 
0.36716E-04 
0.35840E-04 
0.30563E--04 
0.21624E--04 
0. 17163E-04 
I  =  46 
0.46329E-06 
0.45794E-06 
0.38994E-06 
0.26270E-06 
0.  14679E-06 
I  =  49 
0.26242E-07 
0.26203E-07 
0.22283E-07 
0.14626E-07 
0.73123E-08 


0.66251E  00 
0.60629E  00 
0.60448E  00 

0.  50144E  00 
0.48878E-00 
0.41860E-00 
0.  38468E-00 
0. 38349E-00 

0.42898E-00 
0. 42042E-00 
0. 36064E-00 
0.33175E-00 
0.33104E-00 

0.66455E-0 1 
0.647  19E-0  1 
0.54687E-01 
0.49643E-01 
0.49427E-01 

0.20843E-02 
0.20233E-02 
0.  16925E-02 
0.  13732E-02 
0.  12997E-02 

0.36668E-04 
0.35160E-04 
0.29204E-04 
0. 19769E-04 
0.  14458E-04 

0.  4664  6E-06 
0.4491 8E-06 
0.37092E-06 
0.23310E-06 
0.  10272E-06 

0.26572E-07 
0.25704E-07 
0.21 132E-07 
0. 1284 lE-07 
0.57962E-08 


0.50115E  00 
0.47563E-00 
0.40805E-00 
0. 38375E-00 
0. 38375E-00 

0.42896E-00 
0.40936E-00 
0.35150E-00 
0.33121E-00 
0.33 109E-00 

0. 66422E-01 
0. 62880E-01 
0.52798E-0 1 
49463E-01 
3.49484E-01 

0.20874E-02 
0. 19663E-02 
0. 1 6169E-02 
0.  13084E-02 
0. 13385E-02 

0-36591E-04 
0.34  184E-04 
0.27668E-04 
0.  1  7291E-04 
0- 1 61 77E-04 

0.46741E-06 
0. 4371 0E-06 
0.34888E-06 
0.20434E-06 
0.  12604E-06 

0.26712E-07 
0.25026E-07 
0.  19790E-07 
0. 1 1 166E-07 
0.52708E-08 


0.62765E  00 
0.60500E  00 
0.60293E  00 

0.50055E  00 
0.45765E-00 
0.39739E-00 
0.38398E-00 
0.38304E-00 

0.42876E-00 
0.39324E-00 
0.34271E-00 
0.331 18E-00 
0.33092E-00 

0.66345E-0  1 
0.604 17E-0 1 
0.5  1404E-0 1 
0.49523E-0 1 
0.49357E-01 

0.20853E-02 
0.  18958E-02 
0. 15350E-02 
0. 13674E-02 
0.  12247E-02 

0.36465E-04 
0.330 16E-04 
0.25865E-04 
0.  17312E-04 
0. 12024E-04 

0.46634E-06 
0.42265E-06 
0.32288E-06 
0. 18206E-06 
0.75264E-07 

0.26680E-07 
0.24208E-07 
0. 18221E-07 
0.97188E-08 
0.27290E-08 


0.61566E  00 
0.60457E  00 
0.60293E  00 

0.499 15E-00 
0.44301E-00 
0.39018E-00 
0.  38355E-00 
0. 38304E-00 

0.42801E-00 
0. 38093E-00 
0. 33634E-00 
0.33105E-00 
0. 33092E-00 

0.66166E-01 
0.58357E-01 
0.50464E-01 
0.49428E-0 1 
0.49357E-01 

0.20775E-02 
0.  18281E-02 
0.  14389E-02 
0.  13066E~02 
0.12247E-02 

0.36245E-04 
0.31806E-04 
0.23772E-04 
0. 14742E-04 
0.  12024E-04 

0.46342E-06 
0.40691E-06 
0.29344E-06 
0.  15386E-06 
0.75264E-07 

0.26518E-07 
0.23294E-07 
0. 16461E-07 
0. 83054E-08 

0. 
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R  *  0.92437E-01 


UNCLASSIFIED 


RUN  G  4 
TIME  =  0.0028900 
I  =  1 

0.  11769E-05  0.  in92E-05 

0.89472E-06  0.82152E-06 

0.50722E-06  0.43553E-06 

0.21716E-06  0.18019E-06 

0.8964IE-07  0.81537E-07 

I  =  4 

0.50623E-05  0.48912E-05 

0.40289E-05  0.37209E-05 

0.23651E-05  0.20489E-05 

0.10635E-05  0.89145E-06 

0.43985E-06  0.45465E-06 

I  =  7 

0.56059E-04  0.561  79E-04 

0.48656E-04  0.45306E-04 

0.29991E-04  0.26309E-04 

0. 14477E-04  0.12218E-04 

0.82571E-05  0.10701E-04 

I  =  10 

0.64285E-03  0.62797E-03 

0.S3780E-03  0.50377E-03 

0.  34698E-03  0.  30846E-03 

0.19050E-03  0.16726E-03 

0.17229E-03  0.20490E-03 

I  =  13 

0.58879E-02  0.58852E-02 

0.52709E-02  0.49940E-02 

0.36596E-02  0.33312E-02 

0.29829E-02  0.27581E-02 

0.28029E-02  0.30844E-02 

I  =  16 

0.51434E-01  0.50821E-01 

0.46153E~01  0.  44230E-01 

0.36042E-01  0.36481E-01 

0.34719E-01  0.33497E--01 

0.33676E-01  0.35097E-01 

I  *  19 

0.37176E-00  0.37126E-00 

0.35016E-00  0.33952E-00 

0.31929E-00  0.31680E-00 

0.30274E-00  0.29854E-00 

0.29811E-00  0.30195E-00 

I  =  22 

0.22555E  01  0.22484E  01 

0.22105E  01  0.21889E  01 

0.20637E  01  0.20369E  01 

0. 19585E  01  0. 19483E  01 

0.  19423E  01  0.  19464E  01 

I  =  25 

0.60548E  01  0.60728E  01 

0.60588E  01  0.60201E  01 

0.57607E  01  0.57010E  01 

0.55435E  01  0.55300E  01 

0.55172E  01  0.55174E  01 

I  *  28 

0.21137E  01  0.21121E  01 

0-20856E  01  0.20722E  01 


0. 10691E-05  0.10186E-05 

0.74352E-06  0.66340E-06 

0.37014E-06  0.31184E-06 

0.  14947E-06  0.  12438E-06 

0.73029E-07  0.84173E--07 

0.47245E-05  0.45345E-05 

0.33890E-05  0.30447E-05 

0.  17579E-05  0.  14958E-05 

0.74732E-06  0.62469E-06 

0.40320E-06  0.62217E-06 

0.55478E-04  0.53867E-04 

0-41630E-04  0.37764E-04 

0.22877E-04  0.19743E-04 

0.  10720E-04  0.80491E-05 

0.96980E-05  0.13367E-04 

0.61192E-03  0.59225E-03 

0.46645E-03  0.42707E-03 

0.27206E-03  0.23897E-03 

0.19054E-03  0.16690E-03 

0.19289E-03  0.23a02E-03 

0.58230E-02  0.56929E-02 

0.468458-02  0.43523E-02 

0.29047E-02  0.30517E-02 

0.29500E-02  0.27462E-02 

0.29953E-02  0.33592E-02 

0-50077E-01  0.49084E-01 

0.42053E-01  0.39677E-01 

0.34899E-01  0.35429E-01 

0.34384E-01  0.33350E-01 

0.34725E-01  0.36491E-01 

0. 36907E-00  0.36479E-00 

0.33780E-00  0.33014E-00 

0.30948E-00  0.30810E-00 

0-30032E-00  0.29730E-00 

0.301  12E-00  0.30593E-00 

0.22387E  01  0.22234E  01 

0.21641E  01  0.212858  01 

0.20075E  01  0.19878E  01 

0. 19472E  01  0.19424E  01 

0.19453E  01  0.19506E  01 

0.60838E  01  0.60871E  01 

0.59626E  01  0.58917E  01 

0.56473E  01  0.56016E  01 

0.55232E  01  0.55195E  01 

0.55174E  01  0.55182E  01 

0.21054E  01  0.20934E  01 

0.20470E  01  0.20211E  01 

UNCLASSIFIED 


0.96084E-06 
0.58383E-06 
0.26090E-06 
0.  T0483E-06 

0. 

0.43025E-05 
0.26999E-05 
0.  12644E-05 
0.  53695E-06 
0.62217E-06 

0. 51550E-04 
0.33844E-04 
0.  16933E~04 
0.98333E-05 
0. 13367E-04 

0.56753E-03 
0.38685E-03 
0.20512E-03 
0.  19241E-03 
0.23802E-03 

0.55064E-02 

0.40074E-02 

0.28068E-02 

0.29668E-02 

0.33592E-02 

0.47779E-01 

0.37564E-01 

0.34012E-0! 

0.34447E-01 

0.36491E-01 

0.35853E-00 
0.32737E-00 
0. 30236E-00 
0.3001  lE-00 
0.30593E-00 

0.22181E  01 
0.20981E  01 
0.19675E  01 
0.  19447E  01 
0.  19506E  01 

0.60794E  01 
0.58239E  01 
0.55662E  01 
0.55181E  01 
0.55182E  01 

0.20952E  01 
0. 19902E  01 
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RUN  G 

0.  19650E  01 
0.18675E  01 
0.  18547E  01 
I  =  31 
0. 1161  IE  01 
0.11260E  01 
0.10630E  01 
0.10162E  01 
0.  10107E  01 
I  =  34 
0.82654E  00 
0.81452E  00 
0.76234E  00 
0.73338E  00 
0.73020E  00 
I  =  37 
0.22789E-00 
0.222S4E-00 
0.20344E-00 
0.  18967E-00 
0.  1  8874E-00 
I  =  40 
0.35814E-0  1 
0.35089E-01 
0.31509E-0 1 
0.26543E-0  1 
0.26494E-01 
I  *  43 
0.49096E-02 
0.46893E-02 
0.41382E-02 
0.33163E-02 
0.28009E-02 
I  =  46 
0.56278E-03 
0.  552  55E-03 
0.48146E-03 
0.36705E-03 
0.24066E-03 
I  *  49 
0.16086E-03 
0. 16225E-03 
0. 14052E-03 
0.  10348E-03 
0.62368E-04 


4 

0.  19367E  01 
0. 18632E  01 
0.  18510E  01 

0.  1 1  576E  01 
0.11097E  01 
0.  10490E  01 
0.  10152E  01 
0.  10076E  01 

0.82678E  00 
0.80633E  00 
0.75380E  00 
0.73273E  00 
0.72841E  00 

0.22  73  lE-00 
0.2 1986E-00 
0. 19886E-00 
0.  18983E-00 
0.  18760E-00 

0.35893E-01 
0. 34586E-0 1 
0.30574E-01 
0.26786E-01 
0.25947E-01 

0.48710E-02 

0.46109E-02 

0.39878E-02 

0.31406E-02 

0.26548E-02 

0.56697E-03 
0.54283E-03 
0.46109E-03 
0.34  197E-03 
0.20432E-03 

0. 16395E-03 
0. 15943E-03 
0. 13404E-03 
0.95260E-04 
0.53885E-04 


UNCLASS  IF 


0.  19158E  01 
0.18570E  01 
0.  185 17E  0 1 


0.18936E  01 
0.18572E  01 
0.18473E  01 


0.  1  153.3E  01 
0.10985E  01 
0.  10402E  01 
0.  101 1 3E  01 
0.  10082E  0  1 

0.82598E  00 
0.7953bE  00 
0.  74829E  00 
0.73042E  00 
0.72885E  00 

0.22666E-00 
0.2 1637E-00 
0.  1  9445E-00 
0.  18866E-00 
0. 1 8787E-00 

0.35865E-01 
0.33956E-01 
0.29605E-01 
0.26341E-01 
0.2608 lE-01 

0.48363E-02 
0.45155E-02 
0. 38281 E-02 
0.29590E-02 
0.26987E-02 

0.56774E-03 
0.53080E-03 
0.43916E-03 
0.31672E-03 
0.21 122E-03 

0- 16547E-03 
0. 15583E-03 
0. 12699E-03 
0.87018E-04 
0.45917E-04 


0. 1 1471E  01 
0.10882E  01 
0.10207E  01 
0. 10123E  0  1 
0.10042E  01 

0.82384E  00 
0.78228E  00 
0.741 13E  00 
0.73105E  00 
0.72659E  00 

0.22578E-00 
0.21227E-00 
0.  19204E-00 
0.  18920E-00 
0.  18642E-00 

0.357 19E-01 
0.33217E-01 
0.28620E-01 
0.26666E-01 
0.25359E-01 

0.47991E-02 
0.44039E-02 
0.366  16E-02 
0.28231E-02 
0.25031E-02 

0.56  509  E-0  3 
0.51648E-03 
0.41595E-03 
0.29132E-03 
0. 18098E-03 

0. 16540E-03 
0. 15146 E-03 
0. 1 1947E-03 
0.78801E-04 
0.34796E-04 


0.  18805E  01 
0.18535E  01 
0.18473E  01 

0.  1  1383E  01 
0.  10734E  01 
0.  10234E  01 
0.  10094E  01 
0.10042E  01 

0.a2017E  00 
0.76866E  00 
0.73775E  00 
0.72939E  00 
0.72659E  00 

0.22446E-00 
0.20791E-00 
0. 19 145E-00 
0.  18821E-00 
0.  18642E-00 

0.35463E-01 

0.32395E-01 

0.27633E-01 

0.26222E-01 

0.25359E-01 

0.47514E-02 
0.42776E-02 
0.34903E-02 
0-271 64E-02 
0.25031E-02 

0.55997E-03 
0.49999E-03 
0.39180E-03 
0.26545E-03 
0. 18098E-03 

0. 16427E-03 
0. 14634E-03 
0. 1 1  159E-03 
0.70604E-04 

0. 


43 


BIBLIOGRAPHY 


1.  Geeraets,  W.J.,  et  al,  "The  Loss  of  Light  Energy  in  Retina  and 
Choroid”,  Archives  of  Opthamology,  64; 606,  1960. 

2.  Armstrong,  B.H.,  et  al.  Absorption  Coefficients  of  Air  from  22,000° 
to  220,000°,  Headquarters  AFSWC,  Kirtland  AFB,  NMEX,  SWC  TR  58-36, 

15  August  1958. 

3.  Erode,  H.L. ,  and  R.E.  Meyerott,  Thermal  Radiation  from  Atomic 
Detonations  at  Times  Near  Breakaway  (U) ,  The  RAND  Corporation, 

RM-1851,  7  August  1956,  SECRET  RESTRICTED  DATA. 

4.  Erode,  H.L. ,  Close-in  H-Bomb  Effects  (U) ,  The  RAND  Corporation, 
RM-1583-1,  Revised  3  February  1956,  SECRET  RESTRICTED  DATA. 

5.  Brode,  H.L. ,  Point  Source  Explosion  in  Air,  The  RAND  Corporation, 
RM-1824-AEC,  3  December  1956. 

6.  Brode,  H.L. ,  and  F.R.  Gilmore,  Estimates  of  the  Thermal  Radiation 
from  Nuclear  Weapons  Burst  at  High  Altitudes  (U) ,  The  RAND  Corporation, 
RM-1983,  18  September  1957.  SECRET  RESTRICTED  DATA. 

7.  Brode,  H.L. ,  Nuclear  Blast  Wave  Calculations  Including  Thermal 
Losses  (U) ,  The  RAND  Corporation,  RM-2076,  27  December  1957,  SECRET. 

8.  Brode,  H.L.,  Thermal  Radiation  from  Nuclear  Weapons  (U) .  The  HAND 
Corporation,  S-82,  30  April  1958- 

9.  Brode,  H.L.,  Theoretical  Desorption  of  the  Early  Phases  of  the 
Fireball  for  a  Very  High  Altitude  Kiloton  Explosion  (U) ,  The  RAND 
Corporation,  RH-2247,  2  September  1958,  SECRET  FORMERLY  RESTRICTED 
DATA. 

10.  Brode,  H.L. ,  Theoretical  Description  of  the  Blast  and  Fireball  for  a 
Sea  Level  Megaton  Explosion  (U) ,  The  RAND  Corporation,  RM-2248, 

2  September  1959,  SECRET  RESTRICTED  DATA. 

11,  Lerohl,  Major  J.K. ,  Retinal  Burns  from  Nuclear  Detonations  (U), 

The  RAND  Corporation,  RM-2424,  24  July  1959,  SECRET  RESTRICTED  DATA. 

12,  Hillendahl,  R. W. ,  Characteristics  of  the  Thermal  Radiation  from 
Nuclear  Detonations,  Vol  III.  U.S.  Naval  Radiological  Defense 
Laboratory,  AFSWP-902,  30  June  1959,  SECRET  RESTRICTED  DATA. 

13.  Elswick,  W.R.  5  and  H.H.  Mitchell,  M.D. ,  An  Evaluation  of  the  Hman 
Retinal  Burn  Problem  Arising  from  Atomic  Detonations,  The  RAND 
Corporation,  RM-2500,  15  December  1959. 


44 


14.  Meyerott,  R.E, ,  et  al,  Thermal  Radiation  Phenomena  Associated  with 
Nuclear  Explosions  in  Air,  Lockheed  Missiles  and  Space  Division, 

LMSD- 703048,  (DASA  1224),  November  1960,  SECRET  RESTRICTED  DATA. 

15.  Glasstone,  S.,  ed. ,  The  Effects  of  Nuclear  Weapons,  Government 
Printing  Office,  2nd  Edition,  to  be  published  in  April  1962. 

16.  Wray,  J.L. ,  Retinal  Burn  Problem  Evaluation  for  High  Altitude  Nuclear 
Bursts  (U) ,  Headquarters  Defense  Atomic  Support  Agency,  Washington  25, 
D.C. ,  DASA  614,  1  July  1959,  SECRET  FORMERLY  RESTRICTED  DATA. 

17.  Kulper,  G.P. ,  The  Atmosphere  of  the  Earth  and  Planets,  2nd  ed. , 
University  of  Chicago  Press,  Chicago,  1952. 

18.  List,  R.J. ,  compiler,  Smithsonian  Meteorological  Tables,  Smithsonian 
Institution,  Washington,  D.C.,  1951. 

19.  Hoerlin,  H.H. ,  et  al.  High  Altitude  Explosions  and  Eyeburn  Problem  (U) , 
Los  Alamos  Scientific  Laboratory,  LAMS  2337,  August  14,  1959,  SECRET 
RESTRICTED  DATA. 

20.  Informal  communicatlo..  with  Dr.  W.J.  Geeraets,  Medical  College  of 
Virginia. 

21.  Carslwa,  H.S.,  and  J.C.  Jaeger,  Conduction  of  Heat  in  Solids,  2nd  ed. , 
Oxford,  Clarendon  Press,  1959. 

22.  Ham,  W. T. ,  Jr.,  "Flash  Burns  in  the  Rabbit  Retina:  As  a  Means  of 
Evaluating  the  Retinal  Hazard  from  Nuclear  Weapons",  American  Journal 
of  Opthamology,  V.  46,  No.  5,  Pt  1,  November  1958,  pp.  700-723. 

23.  Richtmyer,  F.K. ,  et  al.  Introduction  to  Modern  Physics,  McGraw-Hill 
Book  Company,  Inc.,  New  York,  1955, 


45 


DISTRIBUTION  LIST 


ADDRESSEE 


NO  OF  CYS 


ARMY 

Deputy  Chief  of  Staff  for  Military  Operations,  D/A,  Wash  25,  D.C.  1 

ATTN:  Dir  of  SW&R 

Chief  of  Research  &  Development,  D/A,  Wash  25,  D.C. ,  ATTN:  Atomic  Div  1 

Assistant  Chief  of  Staff,  Intelligence,  D/A,  Wash  25,  D.C.  1 

Chief  Chemical  Officer,  D/A,  Wash  25,  D.C.  2 

Chief  of  Engineers,  D/A,  Wash  25,  D.C. ,  ATTN:  ENGNB  1 

ATTN:  ENGEB  i 

ATTN:  ENGTB  1 

Office,  Chief  of  Ordnance,  D/A,  Wash  25,  D.C. ,  ATTN;  ORDTN  1 

Chief  Signal  Officer,  D/A,  Plans,  Programs  &  Operations  Div.,  Wash  25,  i 
D.C. ,  ATTN;  SIGOP-7A 

Chief  of  Transportation,  D/A,  Office  of  Planning  &  Intelligence,  1 

Wash  25,  D.C. 

The  Surgeon  General,  D/A,  Wash  25,  D.C. ,  ATTN;  MEDNE  1 

Commander- in-Chief ,  U.S.  Army  Europe,  APO  403,  New  York,  N.Y. ,  1 

ATTN;  OPOT  Div.,  Weapons  Branch 

Commanding  General,  U.S.  Continental  Army  Command,  Ft.  Monroe,  Va.  4 

Director  of  Special  Weapons  Development  Office,  Hq  CONARC,  Ft.  Bliss,  1 

Texas,  ATTN:  Captain  Chester  I.  Peterson 

President,  U.S.  Army  Artillery  Board,  Ft.  Sill,  Okla.  1 

President,  U.S.  Army  Infantry  Board,  Ft.  Banning,  Ga.  1 

President,  U.S.  Army  Aviation  Board,  Ft.  Rucker,  Ala.,  ATTN:  ATBG-DG  1 


Commsuidant,  U.S.  Army  Command  &  General  Staff  College,  Ft.  Leavenworth,  1 
Kan.,  ATTN:  Archives 

Commandant,  U.S.  Army  Air  Defense  School,  Ft.  Bliss,  Tex.,  ATTN:  1 

Command  &  Staff  Dept 

Commandant,  U.S.  Array  Armored  School,  Ft.  Knox,  Ky. 

46 


1 


ADDRESSEE 


NO  OF  CYS 


Commandant,  U.S.  Amy  Artillery  &  Missile  School,  Ft.  Sill,  Okla.,  1 

ATTN:  Combat  Development  Department 

Commandant,  U.S.  Amy  Aviation  School,  Ft.  Rucker,  Ala.  1 

Commandant,  U.S.  Army  Infantry  School,  Ft.  Denning,  Ga.,  ATTN:  C.D.S.  1 

The  Superintendent,  U.S.  Military  Academy,  West  Point,  N.Y. ,  1 

ATTN:  Prof  of  Ordnance 


Commandant,  Quartemaster  School,  U.S.  Amy,  Ft.  Lee,  Va.,  ATTN:  Chief  1 
OH  Library 

Commandant,  U.S.  Amy  Ordnance  School,  Aberdeen  Proving  Ground,  Md.  1 

Commandant,  U.S.  Amy  Ordnance  &  Guided  Missile  School,  Redstone  Arsenal  1 
Ala. 

Commanding  General,  Chemical  Corps  Training  Command,  Ft.  McClellan,  Ala.  1 


Commandant,  U.S.  Army  Signal  School,  Ft.  Monmouth,  N.J.  1 

Commandant,  U.S.  Amy  Transportation  School,  Ft.  Eustis,  Va.,  ATTN:  1 

Security  &  Info  Off 

Commanding  General,  The  Engineer  Center,  Ft.  Belvoir,  Va.,  ATTN:  Asst  1 
CMDT,  Engr  School 

Commanding  General,  Amy  Medical  Service  School,  Brooke  Amy  Medical  1 
Center,  Ft.  Sam  Houston,  Texas 


Commanding  Officer,  9th  Hospital  Center,  APO  180,  New  York,  N. Y. ,  ATTN:  1 
CO  U.S.  Amy  Nuclear  Medicine  Research  Detachment,  Europe 

Director,  Armed  Forces  Institute  of  Pathology,  Walter  Reed  Amy  Medical  1 
Center,  Wash  25,  D.C. 

Commanding  Officer,  Amy  Medical  Research  Lab.,  Ft.  Knox,  Ky.  1 

Commandant,  Walter  Reed  Amy  Institute  of  Research,  Walter  Reed  Amy  1 
Medical  Center,  Wash  25,  D.C. 

Commanding  General,  QM  R&D  Command,  QM  R&D  Center,  Natick,  Mass.,  ATTN:  1 
CBR  Liaison  Officer 

Commanding  General,  U.S.  Amy  Chemical  Corps,  Research  &  Development  1 
Command,  Wash  25,  D.C. 

Commanding  Officer,  Chemical  Warfare  Lab.,  Amy  Chemical  Center,  Md. ,  1 

ATTN :  Tech  Library 


47 


ADDRESSEK 


NO  OF  CYS 


Commanding  General,  Engineer  Research  &  Development  Lab,,  Ft,  Bel voir,  1 
Va. ,  ATTNs  Chief,  Tech  Support  Branch 

Director,  Waterways  Experiment  Station,  P.O.  Box  631,  Vicksburg,  Miss.,  1 
ATTN  s  Library 

Commanding  Officer,  Diamond  Ordnance  Fuze  Lab.,  Wash  25,  D.C. ,  ATTN:  1 

Chief,  Nuclear  Vulnerability  Branch  (230) 

Commanding  General,  Aberdeen  Proving  Ground,  Aberdeen  Proving  Ground,  1 
Md.,  ATTN:  Ballistics  Research  Labs,,  ATTN:  Dir,  BRL 

Commander,  Army  Ballistic  Missile  Agency,  Redstone  Arsenal,  Ala.,  ATTN:  1 
OROAB-HT 

Commanding  General,  U.S.  Army  Electronic  Proving  Ground,  Ft,  Huachuca,  1 
Arizona,  ATTN:  Tech  Library 

Commanding  Officer,  USA  Signal  RScD  Lab,,  Ft,  Monmouth,  N.J.,  ATTN:  Tech  1 
Ooc  Center,  Evans  Area 

Director,  Operations  Research  Office,  Johns  Hopkins  University,  6935  1 

Arlington  Rd, ,  Bethesda  14,  Md, 

Commandant,  U.S.  Army  Chemical  Corps  CBR  Weapons  School,  Dugway  Proving  1 
Ground,  Dugway,  Utah 

NAVY 

Chief  of  Naval  Operations,  D/N,  Wash  25,  D.C. ,  ATTN:  0P03EG  1 

ATTN:  OP-75  1 

ATTN:  OP-91  1 

ATTN:  OP-922G1  1 

ATTN;  OP-92262  1 

Chief  of  Naval  Personnel,  D/N,  Wash  25,  D.C.  1 

Chief  of  Naval  Research,  D/N,  Wash  25,  D.C. ,  ATTN:  Code  811  1 

Chief,  Bureau  of  Medicine  &  Surgery,  D/N,  Wash  25,  D.C. ,  ATTN:  Special  1 
Wpns  Def  Div 

Chief,  Bureau  of  Ships,  D/N,  Wash  25,  D.C. ,  ATTN:  Code  423  1 

Chief,  Bureau  of  Yards  &  Docks,  D/N,  Wash  25,  D. C. ,  ATTN:  D-440  1 

Director,  U.S.  Naval  Research  Lab.,  Wash  25,  D.C,  ATTN:  Mrs,  Katherine  1 
H,  Cass 

Commander,  U.S.  Naval  Ordnance  Lab,,  White  Oak,  Silver  Spring,  Md, 

48 


1 


ADDRESSEE 


NO  OF  CYS 


Director,  Material  Laboratory,  New  York  Naval  Shipyard,  Brooklyn  1,  2 

N.Y.  ,  ATTN:  Ml*.  W.  Dirksen  and  Mr.  0.  Mixter,  M.D. 

Commanding  Officer  &  Director,  Navy  Electronic  Lab.,  San  Diego  52,  Calif. 1 

Commanding  Officer,  U.S.  Naval  Mine  Defense  Lab.,  Panama  City,  Fla.  1 

Commanding  Officer,  U.S.  Naval  Radiological  Defense  Lab.,  San  Francisco,  1 
Calif.,  ATTN:  Tech  Info  Div 


Commanding  Officer  &  Director,  U.S.  Naval  Civil  Engineering  Lab.,  Port  1 
Hueneme,  Calif.,  ATTN:  Code  LSI 

Commanding  Officer,  U.S.  Naval  Schools  Command,  U.S.  Navi  ,1  Station,  1 

Treasure  Island,  San  Francisco,  Calif.  / 

i 

Superintendent,  U.S.  Naval  Postgraduate  School,  Monterey,  Calif.  1 

Off leer- in-Charge,  U.S.  Naval  School,  CEC  Officers,  U.S.  Naval  1 

Construction  Battalion  Center,  Port  Hueneme,  Calif. 

Commanding  Officer,  Nuclear  Weapons  Training  Center,  Atljintlc,  U.S.  1 

Naval  Base,  Norfolk  11,  Va.,  ATTN:  Nuclear  Warfare  Dfjpt, 

Commanding  Officer,  Nuclear  Weapons  Training  Center,  Pa‘lfic,  Naval  1 

Station,  San  Diego,  Calif. 

Commanding  Officer,  U.S.  Naval  Damage  Control  Training  C  nter.  Naval  1 
Base,  Philadelphia  12,  Pa.,  ATTN:  ABC  Defense  Course 


Commanding  Officer,  Air  Development  Squadron  5,  VX-5,  China  Lake,  Calif.  1 

Commanding  Officer,  Naval  Air  Material  Center,  Philadelphia  12,  Pa.,  1 
ATTN:  Tech  Data  Branch 

Commanding  Officer,  U.S.  Naval  Air  Development  Center,  .ohnsville.  Pa.,  1 
ATTN:  NAS  Librarian 

Commanding  Officer,  U.S.  Naval  Medical  Research  Institute,  National  1 

Naval  Medical  Center,  Bethesda,  Md. 

Off icer- in-Charge,  U.S.  Naval  Supply  Research  ft  Develops ent  Facility,  1 
Naval  Supply  Center,  Bayonne,  New  Jersey 

Commandant,  U.S.  Marine  Corps,  Wash  25,  D.C. ,  ATTN:  Coce  A03H  1 

Director,  Marine  Corps  Landing  Force,  Development  Center,  MCS,  Quantico,  1 
Va, 


49 


ADDRESSEE 


NO  OF  CYS 


Commanding  Officer,  U.S.  Naval  CIC  School,  U.S.  Naval  Air  Station,  1 

Glynco,  Brunswick,  Ga. 

Chief  of  Naval  Operation,  D/N,  Wash  25,  D.C. ,  ATTN;  0P-09B5  1 

Chief,  Bureau  of  Naval  Weapons,  D/N,  Wash  25,  D.C. ,  ATTN;  RR12  1 

AIR  FORCE 

Hq  USAF  (AFCVC) ,  ATTN;  Operations  Analysis  Office,  Wash  25,  D.C.  1 

Hq  USAF  (AFCIN-3),  Wash  25,  D.C.  1 

Hq  USAFE,  APO  633,  New  York,  N.Y.  1 

Hq  USAF  (AFRDC),  Wash  25,  D.C.  1 

Hq  USAF  (AFCSG),  Wash  25,  D.C.  1 

TAC,  Langley  AFB,  Va.  1 

ADC,  Ent  AFB,  Colo.  1 

AFSC,  Andrews  AFB,  Wash  25,  D.C.  1 

BSD,  AF  Unit  Post  Office,  Los  Angeles  45,  Calif.  1 

PACAF,  APO  953,  San  Francisco,  Calif.  1 

AFCRL,  L.G.  Hanscom  Field,  Bedford,  Mass.  1 

AFSWC,  Klrtland  AFB,  New  Mexico  1 

AUL,  Maxwell  AFB,  Ala.  1 

Tech  Tng  Cen,  Lowry  AFB  30,  Colo.  1 

USAF  Aerospace  Med  Cent,  ATTN;  Lt  Col  Ralph  Allen,  Brooks  AFB,  Texas  1 

1035th  USAF  Fid  Acty  Gp,  Wash  25,  D.C.  1 

ASD,  ATTN;  Aerospace  Medical  Laboratory  (MRVPV) ,  FOR;  Capt.  1 

L.R.  Loper,  Wrlght-Patterson  AFB,  Ohio 

Director,  USAF  Project  RAND,  via;  Air  Force  Liaison  Office,  The  RAND  1 
Corporation,  1700  Main  Street,  Santa  Monica,  Calif. 

RADC,  Grlfflss  AFB,  N.Y.  1 

FTD,  Wrlght-Patterson  AFB,  Ohio  1 

50 


ADDRESSEE 


NO  OF  CYS 


OTHER  AGENCIES 

Director  of  Defense  Research  &  Engineering,  Wash  25,  DoC. ,  ATTN;  Tech  1 
Library 

Director,  Weapons  Systems  Evaluation  Group,  Room  1E880,  the  Pentagon,  1 
Wash  25,  D.Co 

U.S.  Documents  Officer,  Office  of  the  United  States  National  Military  1 
Representative  -  SHAPE,  APO  55,  New  York,  N.Y. 

Chief,  Defense  Atomic  Support  Agency,  Wash  25,  D.C.  10 

Commander,  Field  Command,  DASA,  Sandia  Base,  Albuquerque,  New  Mexico  1 
ATTN:  FCTG 

Commander,  Field  Command,  DASA,  Sandia  Base,  Albuquerque,  New  Mexico  16 

Administrator,  National  Aeronautics  &  Space  Administration,  1520  "H"  1 

St,,  N.W. ,  Wash  25,  D.C. ,  ATTN:  Mr.  R.V.  Rhode 

SAC  {DAWS),  Offutt  AFB,  Nebr  1 

Commandant.  US  Coast  Guard,  1300  E.  Street,  N.W.,  Wash  25,  D.C,,  1 

ATTN:  COIN) 

ASTIA,  Arlington  Hall  Station,  Arlington  12,  Va.  15 

Los  Alamos  Scientific  Lab.,  P.O.  Box  1663,  Los  Alamos,  New  Mexico,  1 

Report  Librarian 

Chief,  Classified  Technical  Library,  Technical  Library,  Technical  1 

Information  Service,  USAEC,  Wash  25,  D.C. ,  ATTN;  Mrs.  Jean  O'Leary 

Medical  College  of  Virginia,  Department  of  Surgery,  12th  &  Broad  Sts.,  2 
Richmond  19,  Va. ,  FOR:  Dr.  William  T.  Ham 

Los  Alamos  Scientific  Lab.,  P.O,  Box  1663,  Los  Alamos,  New  Mexico,  1 

New  Mexico,  FOR;  Dr.  Kerman  Hoerlin  (JIO) 

Mr.  Laurence  S.  Cohan,  OEG  Representative,  c/o  Staff,  Commander,  Anti-  1 
Submarine  Forces,  Atlantic,  U.S.  Atlantic  Fleet,  U. S.  Naval  Base, 

Norfolk  11,  Virginia 


51 


